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signal generated at flame temperatures is quite weak compared to that 
produced at room temperature, our results indicate that it may provide 
an attractive alternative to the widely used vibrational CARS technique. 

In addition to the spatially resolved single point work, we have 

used multipoint CARS to obtain information from many spatially resolved 

3 

volume elements along a cylindrical line (0.1 x 0.1 x 2.0 mm ). For 
this case, a large-angle, phase-matching geometry was used in which the 
dye laser probe beam and one of the pump beams intersect a sheet of 
pump radiation at an angle of 90°. The CARS signal is dispersed by a 
spectrograph and detected by the vidicon in such a way that spectral 
information is obtained in one dimension and spatial information in 
the other. The CARS signal generated in this manner is extremely weak 
since the CARS intensity is proportional to the square of the inter- 
action length, which is 0.1 mm in the large-angle case as compared to 
3 mm in the small-angle geometry. By using an image intensifier with 
the SIT vidicon, we were able to obtain single-shot multipoint spectra 
of N 2 and C> 2 at room temperature. 

We obtained qualitative information on the instantaneous species 
concentration and temperature at 20 spatially resolved volume elements 
(0.1 x 0.1 x 0.1 mm) along a line. Coflowing gases of either different 
species or different temperatures were observed. We used combinations 
of N 2 , 0 2 (which has a rotational spectrum in the same spectral region 
as N 2 ), Ar, and freon (neither of which has a spectrum in this region). 
We also used coflowing hot and cold N 2 » although the results from this 
were not nearly as striking since we found significant mixing of the 
two flows. The boundaries were much more clear-cut in the case of 
different species at the same temperature. Two major difficulties 
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remain in using the multipoint CARS technique to obtain quantitative 
information on temperature and concentration. First, the CARS signal 
must be adequately normalized to take into account random intensity 
fluctuations in the dye laser spectrum and in the spatial profiles of 
the input beams. Second, various distortions in spatial position and 
in intensity are introduced in both the optical path and in the detec- 
tion system (spectrograph, image intensifier, and vidicon) . The 
potential of the multipoint CARS technique to quantitatively determine 
temperature and concentration will depend on the degree to which the 
difficulties with normalization and distortion can be overcome. 



1 . 


INTRODUCTION 


The measurement of temperature and concentration in a gas is an 
important problem in turbulence and combustion research. Accurate 
measurements are required both for diagnostic work in complex systems 
of technological importance (e.g., jet and internal combustion engines) 
and to provide data on simpler systems to aid the development of accurate 
theoretical models. Optical measurement techniques in general are very 
attractive for these applications because they do not require the place- 
ment of a mechanical probe in the medium; they are therefore inherently 
non- intrusive and are potentially applicable in hostile environments 
such as the interior of combustors. Spontaneous Raman scattering in 
particular is recognized as a useful diagnostic technique because Raman 
spectra can be interpreted to yield the temperature and concentration 
of individual components in the scattering medium (Lapp and Penney, 1974). 
Unfortunately, Raman scattering suffers from several drawbacks that 
prevent it from being useful in many applications of interest. The most 
serious difficulties are its low scattering efficiency, which makes 
temporally and spatially resolved measurements difficult, and the fact 
that the Raman-shifted signal is isotropically scattered, which results 
in the weak Raman signals being easily swamped by unwanted background 
light. Coherent anti-Stokes Raman spectroscopy (CARS) is a relatively 
new type of Raman spectroscopy that is not subject to these disadvan- 
tages and at the same time retains many of the best features of spon- 
taneous Raman scattering. 

CARS is a nonlinear optical process that was first discussed and 
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demonstrated not long after the implementation of the first laser in 
1964 (Maker and Terhune, 1965), but has only become a viable diagnostic 
technique in gaseous media in the past several years, due largely to 
the development of reliable high power solid-state and tunable dye 
lasers. Very briefly, CARS is a nonlinear mixing process in which laser 
beams at two different frequencies are combined to produce a Raman signal 
at a new frequency. The CARS process is characterized by high conversion 
efficiency (many orders of magnitude greater than spontaneous Raman 
scattering) and by the coherent, laser-like quality of the output 
signal. These features make CARS applicable to systems in which time- 
resolved measurements are required and to systems with high levels 
of background illumination, since the collimated signal beam can easily 
be spatially filtered to reject the unwanted background light. Like 
spontaneous Raman scattering, the CARS technique generates spectra 
that can be interpreted to yield the temperature and concentration 
of a specific species (Tolies et al., 1977). 

The goal of the work presented in this report was to apply the CARS 
technique to the instantaneous measurement of temperature in nitrogen 
gas with good spatial resolution. The primary temperature range of 
interest was room temperature and below, which is important in wind 
tunnel and shock-tube diagnostic work where large temperature variations 
can be found in regions of rapid expansion. Although CARS is currently 
being widely applied to the measurement of temperature and pressure 
in gaseous systems, most of the work is being done with the large 
frequency shifts associated with the rotational-vibrational Raman lines 
(e.g., the Q-branch in ^ with a shift of 2331 cm *) for application 
to high temperature (1000 to 2500 K) combustion systems (Eckbreth et al.. 
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1979) . At room temperature, however, there is virtually no population 
outside the ground vibrational state, so all temperature information 
must come from the distribution of population among the rotational levels. 
In the present work the rotational populations were deduced from the CARS 
spectra of the pure rotational Raman lines. In principle, the same 
information could be obtained from the rotational lines in the first 
Q-branch, but the pure rotational lines are better resolved and have 
larger Raman scattering cross sections. Extension of the CARS technique 
to the small shifts associated with the rotational Raman lines 
introduces experimental problems not encountered with large shifts, 

so that until recently the only rotational CARS work reported was 

done on the J = 3 -*-5 transition in at a Raman frequency of 1033.4 cm - * 

(Barrett, 1976) and in air down to shifts of 40 cm * (Beattie et al., 

1978) . Very recently, there have been reports of rotational CARS spec- 
tra in the range of 0 to 200 cm * (Goss et al., 1980; Shirley et al., 

1980) , although no extensive attempts to deduce the rotational tempera- 
ture from the spectra have been reported. 

1.1 General CARS Theory 

Coherent anti-Stokes Raman spectroscopy is based on a nonlinear 

optical effect more generally known as four-wave mixing, which is one 

of many nonlinear phenomena that are derived from the third-order, 

(3) 

nonlinear susceptibility, X • Other third-order effects include third- 
harmonic generation, self-focusing and stimulated Raman scattering. All 
of these can be effectively described from a semiclassical viewpoint in 
which optical frequency fields are treated classically with Maxwell's 
equations, and the macroscopic, classical optical properties (e.g., the 
susceptibility) of the nonlinear medium are derived from microscopic 
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calculations utilizing quantum mechanics (Bloembergen, 1965) . This 
section presents a review of the theory of the CARS effect. 

1.1.1 Macroscopic Nonlinear Optics 

The polarization induced in a dielectric medium can be written 
as a power series in the applied electric field as follows: 


- x (1> :E(r,t) + x ( 2 ) :S£.t)E(?,t) 


(l.D 


+ X ' :E(r,t)E(r,t)E(r,t) + ... 


P(r,t) is the macroscopic polarization induced at position r, and 

is a complex, (i + l) th rank tensor. Rather than work in the time 
domain with the nonlinear polarization given in Eq. (1.1), it is cus- 
tomary to work in the frequency domain with the Fourier transform. The 
third-order polarization is then written (Anderson and Hudson, 1978). 


P^ 3) (7,u> 4 ) = D X y' £ ( V , V“ 2 ’“ 3 )Ej (*’“!) 


( 3 ) ,, . 


( 1 . 2 ) 


x E k (r,m 2 )E £ (r,u) 3 )<$(w 4 - - « 2 - Uj) t 


where we have used E(r,t) 


E(r,m) 


-irnt 

e 


du> and have assumed that the 


J 

_oo 

input electric field E(r,u>) is composed of three discrete frequency 
components. Equation (1.2) presents the usual definition of the third— 

(3) 

order, nonlinear susceptibility, X • D is a degeneracy factor that 
takes on the values 6, 3, or 1, depending on whether none, two, or all 
of the frequencies tu^, u> 2> are degenerate. For the type of four- 
wave mixing that is responsible for CARS, the nonlinear polarization 
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+ u >2 - w 3 and is written 


of interest has frequency co. 


p| 3 ) (?,oj 4 ) = 6 x^ £ (oj 4 ; io 1 ,to 2 ,-a) 3 )E j (?,a) 1 )E k (r,a) 2 ) 


( 1 . 3 ) 

* -+ 

X E^(r,to 3 >6(o) 4 - Uj. - ||> 2 + u 3 > . 

(3) 

Because X is a fourth-rank tensor, it must in general contain 

81 elements. In isotropic media, it can be shown that contains 

only 4 unique nonzero elements, X^212’ x 1221’ x 1122’ 3 of whlch 

( 3 ) ( 3 ) ( 3 ) /"a'\ 

are independent since x^ii = X 1212 + x 1221 + x 1122‘ The two Polarization 
cases usually of interest in practical CARS applications are: 


Case (1) 


(1) E(ai 1 ) | |E(ai 2 ) I |E(u; 3 ) , in which 

(’ 3 *) ( 3 ) * 

p ± (w 4 ) = “ 1 > w 2 , ” a, 3^ E i^ a) l^ E i^ a) 2^ E i^ t0 3^ 


( 1 . 4 ) 


Case (2) (2) li (u^) j | E(w 2 ) { E^), in which 

/Q) / 3 } ^ 

P i (“ 4 ) = 6x 1221^“4 ; w i» U) 2’ -W 3^ E j ^l^ E j ( w 2^ E i^3^ * (l- 5 ) 

Although other polarization conditions in which EC^) and E(w 2 ) are 
not parallel are sometimes employed to increase the CARS detectivity 
(Oudar et al., 1979; Rahn et al., 1979), the only polarization condition 
used in the present work is that given in case (2). 

If the electric fields are assumed to be plane waves propaga- 
ting along the z-axis , the Maxwell's wave equation can be reduced 
to the form 


9E (r,a> 4 ) 
3z 


i 


2rrk, 

4 



?(r,u) 4 ) , 


( 1 . 6 ) 
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where we have used E(r,t) = E (r,m)e ^ z ) anc j have assumed that 

"+ * “f ^ 

E (r.oj) is a slowly varying function of r (constant over distances 
^A). Assuming undepleted (constant) input fields E^), E(u 2 > and 
E(w 3 ), Eq. (1.6) can be integrated directly using Eq. (1.5) to obtain 


2-nk, 


/ n \ 

E i( r »^) = 2 ~ f 6 x 1221^4 ; w i * “ 2 ’ ”“ 3 ^ E j (*»“]_)£ j(r,u> 2 ) 


* ^ i-5kz 

E i (r '“3> < >• 


(1.7) 


where Ak - k^ - (k^ + k 2 - k^) . A more general treatment in which the 
fields are not assumed to be copropagating yields the usual vector 
phase-matching relation Ak = - (^ + ^ - £ ) . The detected 

intensity 1 ( 10 ^) is therefore given by (using I = nc/8r J E | 2 ) 


I (“4) 


(8tt)' 


n l n 2 n 3 n 4 



6x 


(3) 

1221 


2 

I(u 1 )I(u 2 )I(a) 3 ) 


x 2 r sin(Ak£/2) , 2 
£ 1 Ak£/2 J ’ 


( 1 . 8 ) 


where Ak£ = Ak*£ and £ is the vector defining the path of signal 
propagation through the nonlinear medium. The integral along the 
direction £j j£^ ^assumed to be along the z-axis in deriving Eq. (1.7)] 
is simply interpreted as the summation of the fields radiated by the 
polarizations at different sites. When the summation is performed in 
the phase-matched direction the fields constructively interfere, 
yielding a large, highly directional signal. For perfect phase matching, 

2 

[ sin(Ak£/2) / (Ak£/2) ] = 1, and the signal intensity increases quadrati- 

cally with pathlength through the medium. The important features 
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exhibited by Eq. (1.8) are summarized "here: 


“4 = “l + “2 ' “3 






I(oj^) « I(cJ 1 )I(a) 2 )I(aj 3 ) 



I(a> 4 ) « l 2 (Ak£ = 0). 

1.1.2 Nonlinear Susceptibility Due to Raman-Type Resonances 

The nonlinear susceptibility can be written as a sum of resonant 
and nonresonant terms. The resonant terms are responsible for the 
Raman-like features in a CARS spectrum and are discussed more fully 
below. The nonresonant term arises from the mixing of the optical 
fields due to the nonlinear response of the electrons in the medium. 

The nonresonant signal generated at is spectrally featureless and 
is usually small compared to the resonantly enhanced signals. An 
energy level diagram that illustrates the mixing process is shown in 
Fig. 1.1(a); note that none of the indicated transitions terminate in 
real excited states and that the initial and final states of the medium 
are identical. This diagram clearly indicates that the four-wave mixing 
responsible for CARS is a parametric process in which energy is trans- 
ferred among the optical fields but is not added to or removed from 

(3) 

the medium. The nonre sonant susceptibility (X^ ) Is, therefore, a 
real quantity, since an imaginary component would imply some loss or gain. 
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Resonant enhancement of the susceptibility can occur whenever one 
of the frequencies to^, ^3 or an Y combination of them corresponds 
to a real transition in the medium. The CARS process results when a 
difference frequency corresponds to a Raman transition in the material. 
The appropriate energy level diagram is shown in Fig. 1.1(b), where 
= has been assumed for definiteness. 

Following the general arguments of the Placzek polarizability 
approximation (Weber, 1973), the electronic polarizability of an indi- 
vidual molecule (the following development implicity assumes that the 
nonlinear medium is a gas) is written (Giordmaine and Kaiser, 1966) 



OL . . 


(a ).. 
o 13 


<t^> 


ij 


where (a ) . . is the polarizability of a molecule with fixed nuclei 
o xj 

and ( 9 a/ 9 q)._. is the change in electronic polarizability due to a 
change in the nuclear conformation q (q can be interpreted as the 
amplitude of the nuclear motion corresponding to a rotational- 
vibrational transition). Thus, the microscopic induced polarization 
in the presence of an applied field ^ is 

p = ( f£) <q> E , (1.9) 

d q 

where a , which leads to elastic scattering, has been neglected. By 
o 

a simple, semiclassical calculation for a two- level system in which q 
is assumed to be driven at frequency ~ by the applied fields at 

u >2 and oj^, it can be shown that (Maier, 1976) 
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NON-RESONANT 


RAMAN RESONANT 


Fig. 1. 


(a ) 


(b) 


1. Energy level diagrams for nonresonant (a) and Raman- 

resonant (b) four-wave mixing. Real and virtual states 
are indicated by solid and dashed lines, respectively. 





:q> 


< f = > 4 < ^ >„ &.*>?*.* * 

2 2 
0) - (oj_ - 03 .) - ir(o)_ - (D,) 

rov 2 3 2 3 


( 1 . 10 ) 


The meanings of the various terms are as follows: 

a) = transition frequency between the upper and lower 
rov 

rotational-vibrational states. 

A = difference in the occupational probability between the 
lower and upper levels (pop. lower - pop. upper) 

2 

q^ = square of matrix element of operator q coupling the 
two states. 

T = damping term; corresponds to the inverse of the de- 
phasing time (l/T^) in the Bloch equations. 


The 


induced polarization at oj^ = + (a^ - w^) °btained by using 


in Eq. (1.9) to obtain 


,-y . , 9 ct . _1 v 

V r ’ t} = ( 9^ > << 1 > E m (r » t) ■ 

4 £m 


(1.11) 


, i(k.z - (o.t) 

-vi ->- ' v j 2 

Recalling that the electric fields are E (r,t) = E (r, to^e , 

and using Eq. (1.10) for <q>, yields 


N ( 




. > W 8? 8^ \ 4E.(r,a, 2 )E (r, m3 )E m (r,» 4 ) 

n xi £m J 


JJL 


“rov " ( “2 " “ 3 )2 " ir( “2 ' “3 } 


( 1 . 12 ) 


X e 


i(k 1 + k 2 - k 3 - k^)z 


Eq. (1.12) gives the macroscopic polarization and was obtained from 
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Eq. (1.11) by multiplying by the molecular density N; in a solid or 

liquid it would also be necessary to include local field corrections. 

Comparison of Eq. (1.12) with the nonlinear polarization derived using 
(3) 

the x formalism [Eq. (1.3)] yields an expression for the Raman- 

(3) 

resonance enhanced susceptibility. Recalling that x (o)^,oo^,oo 2 ,~io 3 ) = 
(3) (3) 

X NR + X RAM (< V “i» w 2 »"“3*» We find that 


0) 


(x RAM (< V aJ l , ‘V tJ 3' )) ijk£ 


N( 


rov .2,3a. ,3a. 

)q 0i ( H ^ > 

ij k£ 


6^i 


“rov - ( “2 - “ 3 )2 “ ir( “2 " “3 ) 


(1.13) 


Equation (1.11) can be reconciled with the corresponding expression 
obtained in a quantum mechanical treatment by replacing some of the 
classical quantities with their quantum mechanical analogues. First of 
all, the perturbation theory calculation yields a true Lorentzian 
lineshape function (Nibler and Knighten, 1979) which differs slightly 
from that in Eq. (1.13); in the remainder of this work the substitution 


2 

oo 

rov 


2co 

rov 



- ir(co 2 - oo 3 ) 


1 


00 

rov 


(“2 ~ " iF 


(1.14) 


will be used. Here, F is recognized as the half-width at half- 
maximum of the Raman line. Also, the polarizability matrix elements 
(a^ ) ^ between the initial and final levels coupled by the Raman 
transition are easily identified with the polarizability derivative 


^im^if 


<i a„ f> 
1 £m' 


<if> 


q if 


(1.15) 
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The population difference factor. A, that appears in Eq. (1.13) 
indicates that the intensity of a CARS peak depends on the difference 
in population between the initial and final rovibrational states. 

This is quite different from ordinary Raman scattering where the 
intensity is proportional to the initial state population only; the 
difference between the two types of spectra is particularly pronounced 
in the case of pure rotational scattering. This leads to a possible 
source of error in CARS spectra because A may change during the course 

of a laser pulse (Tolies et al., 1977). Away from exact resonance 

(3) 

(cc >2 - 7* w rov^ t * ie non lf near susceptibility is composed of and 

(3) (3) 

the off-resonant tail of the Lorentzian line due to Thus, X 

is almost entirely real and no real transitions in the medium can occur. 

In this case, the CARS process is purely parametric [see Fig. 1.1(a)] and 

(3) 

the population difference associated with Xj^j is constant. As to 2 ~ ^ 
is tuned to resonance, however, there is a possibility of real transi- 
tions between the upper and lower rovibrational states. Although the 
CARS process is still strictly parametric, there are other nonlinear 
processes (e.g., stimulated Raman scattering) that occur simultaneously 
with the CARS and can result in the transfer of population from the 
lower to upper state, thus decreasing A. This saturation effect causes 
a loss of signal intensity and distorts A away from the initial value 
determined by the population distribution at thermal equilibrium. This 
must be avoided if an accurate temperature is to be inferred from a 
CARS spectrum. 

In the CARS process the rovibrational transitions are driven by 
the coherent fields E(uj 2 ) and E(u> 3 ), and the values of <q> at different 
molecular sites are related in phase by the factor exp[i(k 2 - k 3 )z]. 

As before, the detected intensity is given by the square of the sum 
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of the scattered fields from all sites 


I(w 4 ) ° 


v v , ±[(t 4 ' *1 " *1 + t 3 )4 - V ] ,2 

l E 1 ( r »‘V e I 


If the summation is performed along the phase-matched direction, the 
phase factor is zero, the scattered fields interfere constructively 
and the cross terms do not vanish. This results in a large, highly 
directional signal, since for any other direction the fields experi- 
ence destructive interference which prevents the build-up of a large 
amplitude scattered field. The detected signal in the phase-matched 
direction is therefore given by 


I( V H l E(oo 1 )E(co 2 )E*( C o 3 > 1 2 


I <xSi> I' I ( W 1 )I( W 2 )I(‘°3 ) 


(1.16) 


Using Eqs. (1.13) and (1.14) in Eq. (1.16) yields 



N( 


12^i 


) ( 


9o_ 

3q 


, 

aq 


k£ 


‘01 


rov 


- (to 2 - a> 3 ) - i T 


(1.17) 


(3) 

Since depends on the square of the polarizability derivative, it 

can be related to the ordinary Raman scattering cross section. 
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and the detected CARS intensity is seen to depend upon the square of 
the ordinary Raman scattering cross section. 

In addition to slimming over all molecular sites in Eq. (1.17) it 
is also necessary to include the nonresonant susceptibility and resonant 
terms from neighboring Raman lines in the sum. Because the detected 
signal depends on the square of the sum of the scattered fields, it is 
possible to have interferences between nearby resonances and between 
the resonant and nonresonant contributions. Therefore, the net nonlinear 
susceptibility is given by 


(3) 


= X 


(3) 

NR 


+ 1 
i 


r — ( — l A 

~ ( 7T ' ( 77T ) A , 

12k “s j J 

0) - (w_ - (D ) - i r. 

i 2 . 3 l 


(1.18) 


where the sum is now over all Raman resonances. Interference between 
the various terms of Eq. (1.18) is a major source of difficulty with 
CARS and will be discussed more fully for pure rotational spectra in 
Section 4. 
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1.2 Nonlinear Susceptibility due to Rotational Raman Resonances In 
Nitrogen Gas 

The general expression for the resonant nonlinear susceptibility 

developed in the previous section will now be specialized to the case 

of pure rotational Raman transitions in nitrogen gas. The intensity 

factors for the Raman transitions are derived from the polarizability 

matrix elements introduced in Eq. (1.15). The exact form of the 

temperature dependent population difference. A, is derived, including 

the nuclear spin degeneracy factor. The final expression obtained for 
(3) 

^RAM w ^'*' directly applicable to the interpretation of the experi- 
mental, temperature dependent CARS spectra generated in this work. 


1.2.1 Evaluation of the Raman Matrix Elements in N r 


The contribution to the nonlinear susceptibility due to the 
rotational Raman transition J -*■ J + 2 can be written using Eq. (1.19) 


(3) RAM 
vx ijkr J 


, N . . , 9 a . 


, 9a N 2 

. . ( 3q q 01 

n k& 


u j " (u) 2 - V " 1 F J 


where oo^ = [E ( J 4- 2) - E(J)]/fi, A^ = population of level J minus 
population of level J + 2, and r is the HWHM of the J •+■ J + 2 Raman 

J 

line. 


A given J rotational level is (2J + l)-fold degenerate in the mag- 
netic quantum number M (assuming no applied magnetic field). Therefore, 
in order to obtain the susceptibility for the entire J to J + 2 Raman 
transition, it is necessary to sum over all degenerate M levels, leading 
to the substitution: 


NAj | <vJM| | v’ J’M’> | ^ 


(2J 


I l^lr. |v'J'M'>| 2 4^ M ’ . 
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where the initial and final states are specified by their vibrational 
(v) and rotational (J,M) quantum numbers. The population factor is 
written to explicitly show that it is associated with the individual 
matrix elements; since all M states are degenerate, however, it can be 
pulled out of the summation. The sums of these matrix elements have 
been evaluated (James and Klemperer, 1959) using the fact that the 
vibrational and rotational contributions to the polarizability can be 
separated for diatomic molecules, as follows: 



The results of the calculation for S-branch scattering, including the 
orientational average, are: 

l |<vJM|a | v' J'M'> | 2 = -? - K(J + 1) |<vJ|y|v’J '>| 2 
MM' 11 

l |<vJM|a Jv'J'M'>| 2 = jr K(J + 1) | <vJ j y | v' J’ > | 2 , 

MM' 

where K(J) = J(J + 1)/(2J + 1) and y is the anisotropy of the molecular 
polarizability. As expected, the isotropic part of the polarizability 
tensor does not play a role in pure rotational scattering, and the 
pure rotational lines are therefore completely depolarized (the 
isotropic part of the polarizability tensor only contributes to AJ = 0 
transitions). The vibrational matrix elements do show a mild J depend- 
ence, but for the low J values used in this work these terms can be 
neglected. 
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1.2.2 The Population Difference Factor 


The population of a rotational state in thermal equilibrium is 
determined by Boltzmann statistics and by nuclear spin statistics. 
Specifically, the population of a particular state characterized by 
quantum numbers J, M is: 

-E(J)hc/KT 

P(JM) = gj (1.20) 


where 


Therefore, 


g = nuclear spin weighting factor 

Q = rotational partition function 

= l g (2J + 1) e -E ( J > hc / kT 
J=0 

J 'M' 

the population difference factor A^ 

A^ + 2* M ' = A = P(J, M) - P(J + 2, M f 

J ,M J 


is 


> 


= g. 


-E(J)hc/KT 

e (1 


-to hc/KT 

e ) 


where A^ is now written to show that it is independent of M. For the 
N 2 molecule in which states with even J are (+) , s, and states with 
odd J are (-), a, we find that the above argument requires that 
gj * 2 for J even and gj = 1 for J odd. 

The final expressions for the nonlinear susceptibility elements 
used in this work are as follows : 
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4N 


(3) .RAM 360k vJ * V J 


, - 0 ) hc/KT 

, 7 ,| S(J) P(J) (1 - e J ) 


( Y w/ ) 

u llll' 


“j ' (a) 2 ‘ “9 " 1F J 


(X 


(3) ) 

1221 J 


RAM 

J 


9 -co hc/KT 

> i Y vJ v'J'I S(J) P(J) (1 ‘ e > 

360^1 VJ,V J 

COj - (co 2 - « 3 ) - 


where S(J) = (J + 1) (J + 2)/(2J + 1) (2J + 3) and P(J) is as given 
in Eq. (1.20). These expressions can be rewritten in terms of the 
ordinary Raman scattering cross section for the J t ^ 1 Raman line as 


(X 


(3) 

1111 


) 


RAM 

J 


_N 

i2h 



s 


2 


a -u ,hc/kT 

( ^ ) P(J) (1 - e ) 


“j - (a) 2 ' <9 ‘ iF J 


(3) RAM 
u 122l' 

J 


3N 

48fc 



s 


2 


„ -oo hc/kT 

( ) P(J) (1 - e J ) 


U)j - (w 2 - u> 3 ) - ir a 
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2 . 


THREE-DIMENSIONAL PHASE MATCHING 


The phase-matching requirement derived in Section 1.1.1 is 


where 



k l + k 2 




n . 

1 


: 2-irn.A * n. to 
1 1 c 

= refractive index at X 


i * 


N^ has no strong absorption lines in the visible region of the spectrum, 
so n is nearly constant and can be ignored. The simplest phase-matching 
geometry, then, is one in which the three input beams and the signal 
beam are colinear, since | k^ | + |k^| = |k^| + | k^ | . This geometry is 
perfectly phase-matched for all Raman shifts and is easily implemented 
experimentally, but it suffers from two serious drawbacks. First, it 
offers inherently poor spatial resolution since the CARS signal is 
generated over the entire region of overlap of the beams. Although the 
colinear beams are generally brought to a focus, at which point the 
signal conversion efficiency is the highest, significant mixing can still 
occur over a region several centimeters in length. Second, it can be 
difficult to separate the relatively weak signal at from the high power 
input beams at frequencies and In rotational-vibrational 

CARS, the shifts are large and it is possible to separate the beams with 
dichroic beamsplitters, but this approach is impossible in pure rotational 
CARS where the Raman shifts are very small (> 10 cm ^) . 
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2.1 Calculation of Phase-Matching Angle 


There is no general requirement for colinear phase-matching and 
several years ago a scheme was proposed in which the three input beams 
are slightly separated in angle (Eckbreth, 1978) . Because signal is 
generated only from the region in which the beams overlap, this approach 
offers greatly improved spatial resolution over the colinear case. This 
geometry is shown in Fig. 2.1 along with the appropriate phase-matching 
equations. If a particular value of a is chosen, the equations can be 
solved for 6 and e: 


cos6 


k cos a 

L 

k cosa 

u 


e 


sine 


k sinS 
o 


where we have used the usual assignments k^ = k 
k^ = k A - The important point to note from these equations is that for 
small shifts (0), the angular separation of the beams is insignificant, 
since k 'v k 'v. k and 0^0 imply that cos6 ^ cosa and sine ^ sinfi. 

L u A 

As a numerical example, choose 0 = 50 cm ^ and a = 5°; then 6 = 5.01° and 
e = 4.99° assuming a pump laser wavelength of 532 nm. 

Because angular separation of the beams is required in pure rotational 
CARS, it is necessary to use a three-dimensional phase-matching geometry 
in which the dye laser beam is directed to be at an angle out of the plane 
formed by the two pump beams (Prior, 1980; Shirley et al., 1980). The 
geometry and phase-matching equations are shown in Fig. 2.2. Because there 
are 5 angles and only 3 equations, it is necessary to choose two angles and 


= k L* 


k~ = k. 


and 
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x : k^cosa + k^ 


cosa + k_ cosa - k„cosS = 


y: k^sina - k^sina - k g sin6 = 


k . cose 
A 


-k. sine 
A 


Fig. 2.1. The BOXCARS phase-matching diagram and equations 



Fig. 2. 


x: k^cosa + k^cosa - 


y: 


k^sina - k^sina - 


z: kgSinB = k^iny 


2. Three-dimensional phase-: 
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kgCosficosB = k^cosecosy 
kgSinScosg = -k^sinecosy 


diagrams and equations. 


solve for the remaining three. Choosing a and g yields 


cosfi 


k^ cos a - 6 
k cosacosg 

b 


siny 


k sing 

b 


k singcosS 

sine = 

k A cosy 


A particularly simple geometry results when 6 = 0 , in which case one finds 


cosg = 


k cos a - 6 

k cosa 
S 


siny = — ^7 


k sing 


e = 0 


This geometry is shown in Fig. 2.3 and corresponds to that used in the 
present work, with a chosen to be 1.5°, so 0 * 1.503“ and y * 1.496“ 
assuming a shift 0 * 50 cm * . 

2.2 Calculation of Coherence Length 

A potential source of error that arises in this geometry and not in 
a colinear geometry is the possibility of phase mismatch over the 
bandwidth of the dye laser. The values of g and y calculated above were 
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Fig. 2.3. The three-dimensional phase-matching geometry used in 
the experiment. 
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based on a Raman shift of 50 cm while the desired broadband CARS 
spectrum may cover the range from 0 to 100 cm Calculation of the 
phase mismatch as a function of the mistunlng away from a shift of 50 cm 
yields 


~ ^^d ' y~~ = 5 '^ cos ^ + Y) - 1) 

where p is the difference between the actual Raman shift and the value 

of 0 for which g and y were calculated. For a maximum mistuning of 

p = ± 50 cm and the values of g and y calculated above, (Ak*£) = 

± £(.068), which yields a coherence length of £ =46 cm (£ , is defined 

coh coh 

such that Ak£ = tt) . Clearly, for an expected interaction length of 
coh 

only a few millimeters, there will be no significant loss of signal 
due to phase mismatch. If a were chosen to be quite large for greater 
spatial resolution, however, the phase mismatch could become a significant 
source of error. Choosing a = 45° and 0 = 50 cm 1 yields g = 45.2° 
and y = 44.8°, in which case the coherence length for a 100 cm” 1 
bandwidth is only £ coh * 628 pan. 
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3. 


EXPERIMENTAL IMPLEMENTATION OF BROADBAND CARS 


A CARS spectrum contains useful information in both the relative 
amplitudes of spectral features and in their absolute intensities. The 
intensity distribution among the Raman features in a given spectrum 
indicate the temperature of the scattering medium, while the overall 
signal intensity is a measure of the density. Although absolute density 
measurements are impractical, the ratio of signal intensities from 
different media can be used to determine their relative density. Unfor- 
tunately, It is difficult to produce experimental spectra that are free 
of spectral distortions and overall amplitude errors. It is important 
to record and analyze the experimental spectra in such a way as to cor- 
rect these errors, if CARS is to be useful as a quantitative diagnostic 
technique. This chapter is concerned with the problem of obtaining 
quantitatively correct spectral information from experimental CARS spec- 
tra; it contains a discussion of various methods for correcting errors 
in the spectra as well as a detailed description of the experimental 
apparatus used in this work. 

3.1 Broadband CARS with Multichannel Detection 

In a typical CARS experiment, the fields at frequencies and 
are provided by a single pump laser beam, the field at frequency is 
provided by a dye laser beam shifted in frequency by 0 to the Stokes side 
of the pump frequency, and the signal is generated at frequency shift 0 to 
the anti-Stokes side of the pump frequency. In this case the following 
frequency assignments are used: 
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The most commonly used method for generating a CARS spectrum is to hold the 
pump laser frequency (to ) constant and to slowly tune the frequency of the 
narrow-bandwidth, Stokes-shifted dye laser (ujg) over the spectral range 
of interest, thus mapping out the CARS spectrum point by point. Because the 
CARS signal is produced only at frequency u> +6, no monochromator is 
required to disperse the scattered light which can be directly detected 
with a photomultiplier tube. In this case the spectral resolution is 
determined by the bandwidths of the pump and Stokes-shifted lasers, which 
can be made very narrow with relatively little power loss. Spectra have 
been obtained with a resolution of 0.012 cm * (Boquillon and Breger, 1979), 
which greatly exceeds that achievable in conventional Raman spectroscopy. 

The single-channel scanning technique is not applicable to making 
time-resolved measurements in dynamic systems, because several minutes may 
be required to scan the entire spectral range of interest. Accurate 
time-resolved measurements require that a complete CARS spectrum be generated 
in a time that is short compared to that required for the system character- 
istics (temperature and/or density) to change. Although these times may 
be as long as many microseconds, most pulsed lasers operate at maximum 
repetition rates of 10' s of pulses per second, which means that the entire 
spectrum must be generated by a single laser pulse ('vlO ns duration for 
the laser used in this work) . 

A broad, single-shot CARS spectrum can be generated simply by deriving 
the Stokes-shifted beam from a dye laser which has a frequency bandwidth 
comparable to the desired spectral interval to be covered (Roh, et al., 

1976) . In this case, a particular spectral component of the dye output 


at frequency (jJ_ * - 0^ will be directly mapped into the anti-Stokes 


spectrum according to the usual mixing rule = 2ou - co, 


“l ■ “s " “l + e i- 
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Thus, for a monochromatic pump laser there is a one-to-one correspondence 
between the dye laser and CARS spectra. The broad-bandwidth CARS signal 
generated in this way can then be dispersed by a spectrograph and 
detected with an optical multichannel analyzer (OMA) . The spectral 
resolution in this case will be determined by the dispersion of the 
spectrograph, by the signal spot size at the spectrograph entrance slit, 
and by the spatial resolution of the OMA. For a nonmonochromatic pump 
laser, however, the ultimate resolution of the CARS spectrum will be 
limited by the pump laser linewidth. 

The most important source of error in spectra produced by the single- 
channel scanning method is the shot-to-shot fluctuation in the intensity 
of the CARS signal, due to pulse-to-pulse variations in the input laser 
intensities. Signal fluctuations as large as ± 30% are common, so that the 
signal from many laser pulses must be averaged at each dye laser setting 
in order to generate spectra that are free of amplitude errors and spectral 
distortion. Other sources of error that cannot be overcome simply by 
averaging many shots are long term drift in the pump and dye laser intensi- 
ties, and the variation of dye laser efficiency with frequency (Wg) . 

Long term intensity drift affects both the relative and absolute magnitude 
of the spectral features, while the dye laser efficiency variation 
distorts individual spectra without affecting their overall amplitude. 

All of these difficulties can be ameliorated by use of a reference 
normalization scheme, but even the most carefully designed scanning systems 
still require that the signals from several shots be averaged to produce a 
single point in the spectrum. 

The broad-bandwidth, multichannel detection technique is not susceptible 
to spectral distortion due to shot-to-shot intensity fluctuations or to long 
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term drift in the laser intensities since an entire spectrum is recorded in 
a single laser firing. The broad -bandwidth dye laser introduces a new 
source of spectral error, however, due to the sharp, random variations of 
laser intensity with frequency. Figure 3.1(a) shows the spectrum of a single 
dye laser pulse; the sharp structures are due to mode beating among the 
many weakly coupled cavity modes present in the dye laser output. These 
structures are directly mapped into the CARS signal, making quantitative 
interpretation of a CARS spectrum virtually impossible without detailed 
knowledge of the spectrum of the corresponding dye laser pulse. The 
spectra are also distorted by the roll-off of average dye laser power at 
the edges of the bandwidth. Amplitude errors are introduced into the 
multichannel spectra by the same mechanisms as in the single-channel 
technique; i.e., shot-to-shot intensity fluctuations and long term drift. 

For non-time-resolved measurements in which it is possible to average 
the signal from many pulses, multichannel detection is superior to the 
single-channel technique for moderate resolution work. The sharp 
structures in the dye laser spectrum tend to average out after many pulses, 
as shown in Fig. 3.1(b), and since the CARS signal is linear in the dye laser 
intensity, so do the structures in the CARS spectrum. At the same time, 
averaging many pulses reduces the overall amplitude error due to shot-to- 
shot intensity fluctuations, just as is the case with scanned spectra. 
Although both methods of obtaining the spectra are susceptible to errors 
due to the variation of dye laser efficiency with a)g , only the scanned 
spectra exhibit spectral distortions due to laser intensity drift. Finally, 
the multipoint technique is much faster than the single-channel technique 
because of the usual multiplex advantage of being able to generate an 
entire, averaged, broadband spectrum in about the same time required to 
generate a single, averaged point in a scanned spectrum. 
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DYE LASER DYE LASER 

INTENSITY INTENSITY 




Fig. 3.1. Dye laser spectra recorded by the OMA for 1 pulse (a), 
and for average of 15 pulses (b) . 
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3.2 Normalization by a Reference 


Although averaging the CARS signals from a large number of successive 
laser pulses corrects many of the spectral and amplitude errors in single- 
and multichannel CARS spectra, it does not account for systematic errors 
caused by laser intensity drift and the dye laser efficiency dependence 
on co s » In addition, averaging is time consuming in the scanning technique 
and is impossible where time-resolved measurements are required. An 
attractive alternative is to normalize the signal from each shot in such a 
way as to make it independent of both the intensities and spectral character- 
istics of the input laser beams. Several different normalization schemes 
will be discussed in this section, and a complete description of the 
normalization done in this work will be presented. 

The most straightforward approach to normalizing a CARS signal is to 
record the pump and dye laser intensities for each shot, and then divide 
the CARS signal by the appropriate product. For single-channel scanning, 
the normalized CARS signal is 


norm . , „ . 

'A <“L + 9) 


V-L + 6) 


W’V-l - 0) 


where I (<o + 0) is the uncorrected, single-shot CARS intensity at shift 0, 

A J_i 

and I (u> ) and I (w - 0) are the intensities of the pump laser and 

J_j i-j , D J-i 

narrow-bandwidth dye laser, respectively. In the multichannel technique 
the entire dye laser spectrum must be recorded in a single pulse, in which 
case the normalized signal is 


_norm , . Q % 
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V-L + V 
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where I^(a)^ + and Ig - 6^) are the intensities of the spectral 
elements corresponding to a Raman shift of 9^ in the broad-band CARS and 
dye laser spectra, respectively. This normalization must be performed on 
the signal from each resolution element of the OMA. The pump laser 
intensity is the same for all 0 and can be neglected if only spectral 
information is required; if accurate intensities are also desired it must 
be included. 

Although this normalization procedure is fairly straightforward, 
experience has shown that it is inadequate when dealing with nonlinear 
interactions. The output of most lasers consists of many modes that 
fluctuate in intensity during the course of a single pulse and are there- 
fore difficult to resolve using conventional techniques. The exact 
temporal, spectral, and spatial relationships between these modes are of 
critical importance in determining the efficiency of a nonlinear mixing 
process, however, so that division by time-averaged (over the duration 
of the pulse) signal intensities fails to account for the observed signal 
intensity fluctuations. 

A better, commonly used approach to the normalization problem is to 
pick off a portion of the input laser beams and mix them in another 
nonlinear medium. Thus, any inter-mode mixing that takes place in the 
sample should be accurately reproduced in the reference as well. If the 
reference medium is chosen to have no resonant features in or near the 
spectral range of interest, the CARS signal will be generated entirely by 
the nonlinear response of the electrons to the applied fields, so that 

x r £ ( “l + 6) I - e) • 
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( 3 ) 

where x^ is the nonresonant susceptibility and Is independent of 0 . An 
atomic gas such as argon makes a good reference medium; although the 
nonresonant susceptibility is small compared to the resonantly enhanced 
Raman susceptibility, adequate signal can be generated by increasing the 
pressure of the gas. With a reference of this type, the normalized signals 
are 


norm 
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for single- channel scanning, and 
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for multichannel detection. 

This type of normalization is easily performed in single-channel 
experiments, since one simply divides the signal generated in the sample 
by that generated in the reference. With multichannel detection, however, 
this is complicated by the necessity of exactly aligning the signal and 
reference spectra so that the spectral component in a given channel of 
the signal spectrum is divided by the same spectral component of the 
reference spectrum. Two types of spatial error in the detected spectra 
occur to make this task nearly impossible: (1) there may be a slight 

registration error between the two spectra, so that one is shifted by a 
few channels with respect to the other; and (2) spatial distortions within 
the spectra occur due to errors introduced by the concave diffraction 
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grating and OMA so that, even if there is exact alignment between the two 
spectra at one place in the field, they are mismatched in frequency else- 
where. In addition to these registration errors, spatial errors in the 
optics and camera also distort the CARS lineshapes. In practice, the 
normalization that is performed is 


norm 

i A 


<“l + V 


ur f <“l + v 1 * 


* 


where [I (to + 0 )], is the signal from the k^ channel of the OMA and 
A L 1 K. 

corresponds to Raman shift 0^. Spatial distortions make it difficult to 
ensure that i = j and, for a dye spectrum with features as sharp as in 
Fig. 3.1(a), errors of only a few channels will be important. 

Because the goal of the present work is to measure the temperature of 
a specific gas rather than to produce rigorously accurate spectra, it is 
possible to use a different type of normalization procedure in which the 
nonresonant reference gas is replaced by the same resonant gas used in the 
sample region. Thus, for identical conditions of temperature and pressure, 
two identical spectra should be produced, one from the sample and one 
from the reference. Instead of performing a channel-by-channel normalization, 
the spectra are simply compared on the basis of the integrated signal 
from corresponding Raman lines. Spatial distortions and registration errors 
do not constitute a problem in this scheme because the individual rotational 
Raman lines in each spectrum are well resolved. The ratio of signals from 
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a particular line will be 


/rJ(w T + e)de 

FjCP.T) = ^ 

/[i a ( “l + s) l refd6 

Any departure of Fj(P,T) from unity should reflect differences in the 
pressure and/or temperature of the sample and reference gases. The 
integrated CARS signal from a particular line is 

U)j+X 

| I][<<*> L + e)de - | ■ 6) l x(3)(t0 j ■ e) l 2 de » 

03j-X 

th 

where u> = frequency of the J rotational Raman line and x >> T_. De- 
J J 

pending on the relative widths of the Raman line and spectral structure 
in the dye laser, I (u - 6) may show strong fluctuations over the range 
of integration or may be relatively constant. By comparing only the 
integrated line strengths of the sample and reference spectra, however, it 
becomes unnecessary to obtain detailed information about the dye spectrum, 
since this is taken care of in the nonlinear mixing process itself. 

The procedure for determining the sample gas temperature is as follows 
The reference gas is held at fixed temperature and pressure, for which 
the theoretical integrated line- intensity distribution is calculated. 

The single shot experimental reference spectrum will, in general, differ 
from the calculated distribution because of the nonuniform dye laser 
spectrum. Therefore, each of the Raman lines in the experimental 
reference spectrum is multiplied by the appropriate factor to bring it 
into agreement with the intensity of the corresponding calculated line. 
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These same correction factors are then applied to the sample spectrum, 
since a given Raman line in both spectra will depend on the same dye 
laser spectral structure. The corrected, integrated line- intensity 
distribution of the sample gas can then be fit to calculated distribu- 
tions to find the temperature. The principal difficulty with this 
approach is the variation of the Raman linewidth with temperature and 
pressure; this will be discussed in Section 4. 

3.3 Optical Arrangement and Alignment 

The pump laser beam at frequency and the broadband, Stokes-shifted 
dye laser beam at frequency a)g are both derived from a high power, pulsed 
Nd:YAG laser. The laser consists of an oscillator and single amplifier 
stage which emit fundamental radiation at 1.06 ym in 15 nsec pulses with 
a maximum energy of 750 mJ per pulse when operating at the design repeti- 
tion rate of 10 Hz. The 1.06 ym radiation is frequency doubled via a 

* 

second-order nonlinear process in a KD P crystal to produce up to 
250 mJ/pulse at 532 nm. This second-harmonic radiation is mixed with the 
residual fundamental in a second nonlinear crystal to produce the third- 
harmonic radiation at 355 nm. The full-power energy distribution in the 
second- and third-harmonic beams is about 160 mJ/pulse and 110 mJ/pulse, 
respectively. 

The second harmonic of the YAG laser is used as the pump beam in 
the CARS experiment, and its linewidth therefore determines the resolution 
of the experiment. For the present work the Nd:YAG oscillator was operated 
without an intra-cavity etalon, so many longitudinal cavity modes could 
oscillate simultaneously; in this case the manufacturer specifies a 
fundamental linewidth of 0.5 cm * FWHM. This parameter was not measured 
in the present work and is likely to vary from pulse to pulse, but a 
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DF 


ODC 


ADC 



Fig. 3.2. Schematic lay-out of the broadband dye laser constructed 

for the experiment. The symbol code is: BS, beamsplitter; 

R, 100% reflector; CL, cylindrical lens; DF, dielectric 
filter; ODC, oscillator dye cell; ADC, amplifier dye cell; 
A, aperture; FR, front reflector of dye laser cavity. 
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value in this range is consistent with the experimental results. 

The third harmonic is used to pump a broad-bandwidth dye laser that 
was constructed for this project. The dye laser consists of a trans- 
versely pumped oxcillator and single-stage, transversely pumped amplifier 
and is shown schematically in Fig. 3.2. Two separate flow systems were 
used to circulate the coumarin 500 dye (in ethanol) through the dye cells 
since the oscillator and amplifier stages require different dye concentra- 
tions. The bandwidth of the emission was unrestricted by any of the usual 
line narrowing devices (prisms, diffraction gratings, etalons, etc.), while 
the center frequency of the band was adjusted by angle tuning a sharp, 
edge-cut, multilayer dielectric filter placed inside the cavity. When 
tuned to a center wavelength near 532 nm, the bandwidth was approximately 
40 l at the half-power points. The pulse duration of the dye laser was 
about 10 nsec and the maximum attainable energy was ^10 mJ/pulse. 

The experimental arrangement is shown in Fig. 3.3. The horizontally 
polarized dye laser beam, denoted <Ci)g>, and the vertically polarized 
pump beams (oj^) are focused by a 30 cm focal length lens into the sample 
region (S) . The input beams, along with the CARS beam generated in the 
sample, are recollimated by the second 30 cm lens. The CARS beam is 
picked off and directed into the spectrograph, while the pump and Stokes- 
shifted beams are refocused into the reference region where the second 
CARS signal is generated. The two horizontally polarized CARS signals 
are directed through polarizers that reject scattered 532 nm light 
and are then focused to pass side by side through the 100 ym wide hori- 
zontal entrance slit of the spectrograph. The two signals are vertically 
dispersed in the spectrograph and are imaged side by side onto the OMA. 

Alignment of the beams to the correct phase-matching angles is 
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Fig. 3.3. Schematic lay-out of the experimental apparatus. The 
symbol code is: T, telescope; BS, beamsplitter; VP, 

vertical-pass polarizer; L, lens; HP, horizontal-pass 
polarizer; HS, horizontal slit; CG, concave grating; 
S, sample gas region; REF, reference gas region. 
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achieved by making sure that the beams are parallel and properly spaced 
before the 30 cm lenses are installed. With the lenses in place, the 
beams are checked for spatial overlap at the foci with a 100 ym diameter 
pinhole. The telescopes, T, are used to bring the dye and pump beams to 
the same diameter OW mm) and to make fine adjustments on the beam 
divergences to ensure that all three beams are brought to focus at the 
point of overlap. Although the alignment is normally done at low power, 
it is adequate to ensure that a signal will appear when the power is 
increased, at which point the signal intensity can be maximized by making 
fine angular adjustments on the beams. When properly aligned, the signals 
from the sample and reference regions can be maximized simultaneously 
even though there is no provision for individual alignment. 

The series configuration for the sample and reference regions was 
chosen over the parallel approach in which portions of the input beams 
are separated from the main beams and focused separately into the 
reference cell to avoid the loss of signal brought about by splitting the 
beams. When a nonresonant reference is used, only a small fraction of 
the input beams need be focused into the reference cell since the signal 
can be boosted by increasing the pressure in the cell; this is undesirable 
with a resonant reference because the Raman linewidth increases linearly 
with pressure, making comparison to the atmospheric pressure sample diffi- 
cult. Splitting the beams evenly between the signal and reference regions 
would mean a loss of signal by a factor of 8 compared to the series 
approach. In principle, there is a danger of the dye laser spectrum 
being perturbed by the mixing in the signal region, resulting in a 
distorted CARS spectrum from the reference, but in practice the beam 
energies are kept low enough so that this is not a significant effect 
(CARS signal intensity is < 10 ^ of the input dye laser intensity) . 
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The spectrograph used for the room and low temperature work consists 
of a single 1 m focal lengthy concave, holographically ruled grating. It 
has a dispersion of 1530 cm m in the focal plane in the spectral region 
of interest. The spectra could be imaged 1:1 onto the OMA or magnified 
by a factor of two, depending on the desired range of coverage. The 
grating in this instrument is not optimized for use as a spectrograph 
and is probably responsible for most of the distortion observed in the 
spectra on the OMA. Fortunately, use of the resonant reference eliminates 
this as a serious source of difficulty. 

All of the experimental spectra were taken with the sample and 
reference regions flooded with pure ^ gas. Although strong signals 
were easily obtained in air, interference from neighboring 0 2 Raman 
lines made quantitative interpretation of the ^ CARS spectra difficult. 
An additional difficulty with an air environment is the presence of dust 
particles that drift into the laser focal volumes and spark, resulting 
in a bright, broad-band signal at the detector. In order to maintain 
dust-free, 100% environments, open-end tubes were placed coaxially 
with the laser beams at the sample and reference mixing regions. Clean 
was then flowed gently into each tube at its center, thus maintaining 
a steady outward flow of from the tube ends. Thermocouples were 
mounted in each tube near the laser beam waists to monitor the gas 
temperature. For the room temperature measurements, the sample and 
reference gas temperatures agreed to within 1 K without any temperature 
control . 

An open-end cell was also used in the low temperature measurements 
to avoid the use of double windows that would have been required with a 
sealed cell. A schematic of the cell is shown in Fig. 3.4(a). ^ gas 
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from a high pressure tank was passed through copper coils immersed in 
L^, through aim insulated plastic tube, and was then blown across 
the cell and out the opposing port. The cell was constructed of copper 
tubing and fittings and was wrapped in foam rubber to alleviate the 
problem of water condensation at low temperature. A 100% ^ atmosphere 
was maintained in the cell by flooding the ends .of the tubes with room 
temperature ^ gas. Without this Nj purge, air was drawn into the tube 
ends and resulted in scattering of the laser beams from frozen water 
vapor at low temperatures. A copper/constantan thermocouple (3 mil wire) 
was mounted in the cell to monitor the temperature. The cell position 
was adjusted so that the thermocouple junction was M. mm above the 
focal point of the beams (closer positioning resulted in spurious readings 
due to the laser beams) . The temperature at the mixing region was adjusted 
by varying the gas flow rate. Although this system does not provide 
any means of active temperature control, the temperature reading from the 
thermocouple was found to stabilize at a given flow rate long enough to 
make a series of measurements. Data in the range of 135 to 250 K were 
obtained with this system. 

The low pressure measurements were made with the 19.5 cm long opti- 
cal cell shown in Fig. 3.4(b). Single-shot spectra could be obtained 
down to a total pressure of about 40 Torr, so the pressure over the entire 
range of interest was conveniently measured with a mercury manometer. The 
high pressure measurements were made with a stainless-steel cell having 
an optical path length of 100 cm. Pressures above 1 atm were measured 
with a single gauge calibrated in 5 psig increments from 0 to 200 psig. 
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3.4 The Optical Multichannel Analyzer 

The optical multichannel detector (OMA) used in this project is 

a digitally controlled television camera (silicon diode array vidicon) 

with a single-stage intensifier. The maximum scanned area is 

2 

1.25 x 1.25 cm and the highest resolution pixel size is 25 ym in 
diameter, although any number of adjacent pixels can be grouped together 
to establish scan formats of lower resolution. The video putput is 
converted to digital form by a 14 bit analog-to-digital converter and 
is passed on to a microcomputer that is used to establish the desired 
scan format as well as to collect and display the data. 

In the present experiment the camera is scanned in two one- 
dimensional tracks of 499 channels each, as shown in Fig. 3.5. Each 
channel consists of a horizontal strip 150 to 200 pixels wide and 1 
pixel tall. The integrated signal from each channel represents a single 
point in the one-dimensional track. This format provides maximum resolu- 
tion in the vertical direction of dispersion and still collects all of 
the signal spread out in the horizontal direction due to astigmatism in 
the concave grating. The 499 x 1 channel arrays from each track are 
stored separately as independent pictures. The scanning rate selected 
for most of this work was 40 ysec/channel, so the time required for one 
complete sweep of both tracks (corresponding to a single frame) is 
39.9 msec, for a rate of 25 frames per second. 

The characteristic of the OMA that is most troublesome for pulsed 
instantaneous measurements is image lag, an inherent problem with many 
types of television tubes whereby the optical signal impressed upon the 
target is not fully erased by a single sweep of the read beam. This 
characteristic gives rise to two problems: signal holdover and nonlinear 
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Fig. 3.5. The two-track scan format used with the OMA to simultaneously 
record the sample gas and reference gas spectra. Each 
track consists of 499 channels spaced on 25 p centers. 
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response (Crowell and Labuda, 1971). Because the signal is not completely 
erased in one or even several sweeps, the signal recovered from any scan 
will contain contributions from several preceding optical pulses. With 
the laser operating at 10 Hz and the OMA scanning at a rate of 25 
frames/s, there will be only 2.5 scans between laser pulses and signifi- 
cant signal hold over does occur. Thus, the signal from any shot does 
not represent a truly instantaneous spectrum. The second problem has to 
do with the fact that the signal is not read off the target linearly; 
it is a common feature of vidicons that the brighter signals tend to 
come off faster than the weak signals, so that after a given number of 
scans there will be an error in the intensity ratio of a strong spectral 
line to a weak line. This shows up dramatically in the CARS spectrum 
of N 2 in which the intensities of alternate peaks are in the ratio of 
4:1. Figure 3.6(a) shows the integrated signal from a single line as 
a function of intensity with the laser operating at 10 Hz. For these 
conditions it should be noted that the signal holdover from previous 
shots tends to offset the inherent nonlinearity so that the overall 
response is fairly linear over an entire decade; without this accidental 
cancellation, the nonlinearity would be considerably worse. 

The solutions to these two problems are easily implemented but 
tend to degrade the overall performance of the system. The problem of 
signal holdover can be eliminated by increasing the time between laser 
pulses until no significant signal overlap occurs. Although the laser 
flash lamps must fire at a steady 10 Hz rate because of thermal lens 
effects in the laser rods, it is possible to gate the Pockels cell so 
that lasing action only occurs for every N th lamp discharge. It was 
found that allowing the laser to pulse only every 1.6 s (N = 16) was 


46 



OMA SIGNAL OMA SIGNAL 
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INCIDENT LIGHT INTENSITY 


OMA signal as a function of incident light intensity. 

In (a), the OMA was scanned once after each laser pulse 


the laser repetition rate was 10 Hz. In (b) , the OMA 
was scanned 10 times after each laser pulse; the laser 
repetition rate was .63 Hz. 




adequate to ensure that all the signal from one pulse was erased before 
the next pulse arrived. The nonlinearity problem was solved by scanning 
the target many times following the arrival of a signal pulse, and 
summing the signal from successive scans to form a single data array for 
each track. It was found that scanning the target 10 - 15 times was 
adequate to ensure that all the signal was read off the target. The 
camera linearity was greatly improved by employing this scan format, as 
seen in Fig. 3.6(b), which shows the integrated detected signal strength 
as a function of incident light intensity and was obtained using 10 
reread frames. 

The price paid for linear, single-shot detection is a slower data 
acquisition rate and noisier signals. The reread frames ensure linear 
response by extracting all the signal from the target but tend to 
accumulate background counts due to the thermally generated dark noise 
in the target and due to electronic noise introduced by the A/D con- 
version process. Although much of this can be subtracted away, there 
is still a noticeable degradation in the signal/noise of the final 
spectra. 
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4. 


CALCULATION OF ROTATIONAL CARS SPECTRA AND INTENSITY DISTRIBUTIONS 


The temperature-fitting procedure used to analyze the CARS 
spectra requires comparison of the distribution of integrated power 
in the peaks of an experimental spectrum with the theoretical 
distributions calculated over a range of temperatures. Therefore, 
the desired calculated quantities are 


A J = 


I 

J 


V“A )d “A 


VV d “A 


I I 


(4.1) 


where I (oj ) is the CARS power as a function of frequency in the 

A n 

vicinity of the Raman line, and the integration extends midway 

to the adjacent peaks on either side of the peak centered at to^. The 
denominator is simply the total CARS power in all lines and is 
included as a normalization factor to allow comparison with experi- 
mental distributions. This chapter is primarily concerned with the 

calculation of the terms 1^. 

A 

4 .1 Modification of the CARS Theory to Include a Pump Laser with 
Finite Linewidth 

In the CARS theory developed In Section 1.1, it was assumed that 
the applied electric field consisted of three discrete frequency 
components, e.g., the outputs of three perfectly monochromatic lasers. 
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In practice, however, lasers have finite spectral linewidths, in which 
case Eq. (1.8) must be modified as follows: 


I(co 4 ) “ 


/ O \ 

da3 1 da) 2 da3 3 | x ^ (cl> 4 ; 


Cl) 
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Ur 


a> 3 ) 


"I 1 (a^)^ (a) 2 )I 3 (a> 3 ) 


(4.2) 

x 6(a>^ - co^ - a> 2 + o) 3 ) , 


where I (u>^) is the power distribution with frequency of laser i. 

A complete discussion of this equation is given in the literature 
(Yuratich, 1979; St. Peters, 1979; St. Peters, 1980), but the basic 
assumption is that the laser outputs consist of densely spaced, perfectly 
monochromatic, independent modes, in which case it is possible to write 
the total CARS intensity at as the sum of intensities due to all 

combinations of co^ , co 2 , and co 3 satisfying the condition co^ = + co 2 - 

This assumption is not fulfilled in real lasers, but it is a useful 
working approximation; an exact treatment that takes into account the 
time-varying frequencies, divergences and degree of coherence of all 
the laser modes is impossible. 

In the ideal broadband CARS experiment the dye laser is assumed 
to have a flat spectrum over the range of interest, so the integration 
over a)^ can be performed in Eq. (4.2) using ^(t^) = I^ and the delta 
function to obtain 


W ‘ *3 


dw-jdo^jx^ w 1 , u> 2 , 0)^ - U) 1 - (i) 2 ) | Z I 1 (ci) 1 )I 2 (cl) 2 ) 


(4.3) 
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For the special case in which a genuinely monochromatic pump laser 
provides the fields at and u) 2 (laser linewidth « Raman linewidths) , 
Eq. (4.3) can be evaluated using 1^(0^) = 1^6(0^ “ W L ) and = 

12*5 ( 0^2 “ to^) to obtain 

Ku) 4 ) « I 1 I 2 I 3 lx (3) (w 4 ; 0> L , 0> L , o> 4 - 2 o> l ) | 2 , 


in which case the integrated CARS power is simply 


T 4 * W3 


( 3 ) 2 

dw 4 !x k ^(“ 4 ; «*> L » w 4 - 2 oj ) | 


(4.4) 


In general, however, the lasers are not monochromatic and the integrated 
CARS power is given by 
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The nonlinear susceptibility at a given wavelength is the sum of 
all resonant and nonresonant terms, so that 
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and 
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This expression explicitly provides for the possibility of simultaneous 

Raman resonances in (= - Wg) and (= co 2 - u)g) . The 

factors Kj contain temperature-dependent and J-dependent terms as well 

(3} 

as all constants associated with expressions for the Kj's were 

derived from Eq. (1.18). In the case of a monochromatic pump laser, 

W 1 = W 2 ~ an< * E< *' (4.5) collapses to a single term. 

In a true four-color CARS experiment (oj^ 4 io 2 4 oj^) it is unlikely 
that both terms in Eq. (4.5) will be on or near resonance simultaneously, 
in which case the off-resonant term can be lumped in with y . in a 
three-color CARS experiment, however, and o j 2 are simply different 
spectral components within the bandwidth of the same pump laser and both 
terms may be near resonance simultaneously, since - co 2 ‘v (I^ = FWHM 
of the pump laser linewidth; 0.5 to 1.0 cm * in this work). Several 
treatments of the effect of finite laser linewidths on CARS spectra have 
appeared in the literature (Yuratich, 1979; St. Peters, 1979; St. Peters, 
1980), but they include only one of the two possible resonance terms, in 
which case the resonant susceptibility is simply 

k j 

*J ( V V V “4 - “l - "2> * - » 2 > - il-j • 


The integrated line strengths corresponding to these two cases are 
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for the double-resonance theory and 
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for the single-resonance case. Since the Xj ' 's in Eq. (4.7) are 
independent of the integration over can be performed and 
Eq. (4.7) can be reduced to the same expression obtained for the case 
of a monochromatic pump laser [Eq. (4.4)], so that 
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Thus, according to the single-resonance theory, the finite pump laser 
linewidth broadens the detected CARS lines but does not affect their 
integrated area. 

The theoretical integrated intensity distributions as well as 
complete rotational CARS spectra have been calculated on a computer 
using both the single-resonance [Eq. (4.8)] and double-resonance 
[Eq. (4.6)] theories. The calculated distributions obtained with the 
two theories are somewhat different, but not enough to significantly 
influence the results of the temperature fits of the experimental data. 
The two types of calculated spectra also show subtle differences, but 
both correspond closely to experimental spectra and it is not possible 
to definitively prefer one theory over the other. Because the present 
experiments are not able to discriminate between the two theories, and 
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because the computational times required for evaluation of the triple 
integrals in Eq. (4.6) are very long, the simpler, monochromatic theory 
was used to calculate the integrated intensity distributions used in 
fitting the data to temperature. 

The CARS signal produced by a given combination of w , , and 

contains contributions from the nonresonant background and from 
all Raman resonances. Because these terms are all summed coherently, 
it is possible to have interferences between on- and off-resonance 
Raman contributions and between the Raman terms and the real, non- 
resonant background. The susceptibility due to an individual Raman 
line can be broken into real and imaginary parts, as follows: 


K J 

“j " x " iF J 


= K n 


(U J ‘ X) 


<“J - X)2 + F J 2 


1F J 

(“J - X)2 + r j 


The real and imaginary lineshapes are plotted in Fig. 4.1. Although 
the imaginary part of the susceptibility falls off rapidly away from 
line center, the tails of the real part die off slowly and can contribute 
to the real, nonresonant background at nearby lines. The effect of this 
real background is to distort the observed CARS lineshapes away from 
the Lorentzian shapes observed in spontaneous Raman scattering. The 
resulting asymmetric lineshapes can be observed in the experimental 
CARS spectra shown in Figs. 4.2(b) and 4.3(b). 
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CARS INTENSITY 



RAMAN SHIFT 


Fig. 4.2. Calculated (a) and experimental (b) CARS spectra of N2 
at 296 K and 1 atm. The experimental spectrum is an 
average of the spectra from 200 consecutive laser firings. 
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Fig. 4.3. Calculated (a) and experimental (b) CARS spectra of N 2 at 
296 K and 15 atm. The experimental spectrum is an average 
of the spectra from 100 consecutive laser firings. 


4.2 Calculation of Quantitatively Correct CARS Spectra 


In order to calculate Intensity distributions that are quanti- 
tatively correct, the values of K^, Tj, and must all be accurately 
known at the desired condition of temperature and pressure. The factors 
K [see Eq. (4.5)] are principally responsible for the temperature depen- 

J 

dence of the rotational CARS spectra, since they include both the number 
density, N, and the population difference, A. The linewidths Tj vary 
with temperature and pressure; they are important parameters because 
they determine the degree of interference between the various suscepti- 
bility terms. The nonresonant susceptibility is constant but its effect 
on a CARS spectrum does depend on the pressure and temperature of the 
gas. 

The rotational Raman lines in ^ are pressure-broadened down to 
below 100 Torr and T therefore scales linearly with the inverse of 
the mean time between molecular collisions. This implies that the 
pressure and temperature dependence of are 


r 


J 


P 


a 


l/T 


1/2 


(4.9) 


In addition to being functions of T and P, the Raman linewidths also 

depend on J. Table 4.1 lists the measured broadening parameters for 

the pure rotational Raman lines of at 300 K (Jammu et_ al. , 1966). 

Since the ^ lines are pressure-broadened at one atmosphere, these 

parameters are simply the linewidths Tj (HWHM) at STP. The J-dependence 

of the T's is also a function of temperature due to the change in the 
J 

rotational state populations with temperature. Unfortunately, this 
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Table 4.1 


Rotational Quantum Number, J 


Broadening Coefficient 


(cm 


-1 


/atm) 


0 

.058 

1 

.055 

2 

.052 

3 

.051 

4 

.048 

5 

.048 

6 

.046 

7 

.042 

8 

.041 

9 

.040 

10 

.038 

11 

.037 

12 

.035 

13 

.034 

14 

.032 

15 

.032 

16 

.031 


These values were not reported; they are obtained by extrapolation the 
data from Jammu et al., (1966). 
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dependence has not been experimentally measured or theoretically 
calculated for the low J levels in cold Nj* [Some experimental and 
theoretical work on this problem has been done recently for Q-branch 
lines in N 2 gas at combustion temperatures (Rahn et al . , 1980; Hall, 

1980)]. Therefore, this dependence was ignored in this work and the 
linewidths at non-STP conditions of temperature and pressure were calcu- 
lated from the values given in Table 4.1 using only Eq. (4.9). 

The remaining quantity required in evaluating Eq. (4.8) is the 
nonresonant susceptibility; in fact, since absolute signal amplitudes 

are not required, it is sufficient to know only the ratio y /K t . 

NR J 

From Eq. (1.18) it can be seen that this is essentially a matter of 
knowing Xj^/ (do/dfl) . Although x^ was measured in ^ several years ago 
(Rado, 1967), no recent measurements are available and no absolute measure- 
ments of the Raman cross sections of the pure rotational lines in N„ have 
been reported. Therefore, the value of x»ro used in this work was chosen 
on the basis of the best visual fit between experimental and calculated 
spectra, judging largely by the asymmetry of the rotational lineshapes 
arising from the interference between x^ and Xj^* The value of 
X NR /(d a /dn) obtained in this way is about 3.5 times larger than that 
obtained using the value of Xj^ reported by Rado (1967) and assuming 
that the pure rotational Raman cross section is about 10 times larger 
than the vibrational Q-branch cross section. This is in fair agreement 
with the finding by other workers that the value of x^ given by Rado (1967) 
is too small by a factor of 2.5 (A. C. Eckbreth, United Technologies 
Research Center, East Hartford, CT, personal communication). 

Spectra calculated using the double-resonance theory for N» at room 
temperature and pressures of 1 atm and 15 atm are shown in Figs. 4.2(a) 
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and 4.3(a), respectively. For comparison, experimental spectra recorded 
at corresponding conditions of temperature and pressure are shown in 
Figs. 4.2(b) and 4.3(b). The experimental spectra are averages of the 
signals from 200 consecutive laser pulses and have been corrected for 
the variation of dye laser efficiency with wavelength. The pump laser 
linewidth used in calculating the theoretical spectra was 1.0 cm * FWHM. 

The sensitivity of the rotational CARS spectra to a small change in 
temperature is shown in Fig. 4.4. The percentage change in the integrated 
intensity of individual rotational lines for a 1 K change in temperature 
at 300 K is plotted as a function of J. For comparison, the integrated 
intensity distribution at 300 K is also shown. The small change in the 
distribution with temperature and the fact that the strongest CARS lines 
are not the most sensitive to temperature changes indicate the difficulty 
of making accurate temperature fits of the experimental spectra. In 
order to minimize the errors introduced by noise in the data, it is 
therefore desirable to include as many rotational CARS peaks as possible 
in the fitting routine in the hope that the noise will average out. 
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Fig. 4.4. The absolute value of the change in integrated signal of 
the individual rotational CARS lines for a 1 K change in 
temperature at 300 K (solid circles). Also shown is the 
integrated signal distribution at 300 K (open circles) . 
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3 DIFFERENCE 



5. 


EXPERIMENTAL RESULTS 


The results of the single-point CARS experiments are presented in 
this section. The results of the experimental investigation of CARS ther- 
mometry for room temperature and cold Nj are presented in Section 5.1. 

The experimental comparison between rotational and vibrational CARS ther- 
mometry for N 2 at flame temperatures is presented in Section 5.2. The 
simultaneous observation of rotational CARS and CSRS for air and room 
temperature N 2 is presented in Section 5.3. 

5. 1 Room Temperature and Cold ^ 

This section presents the experimental results for room temperature 
and cold N 2 - We outline the reduction of the signals recorded on the OMA 
to accurate representations of CARS spectra that can be compared to 
theory. The reduction procedures include several corrections that com- 
pensate for the nonideal characteristics of the detection system. The 
results of determining temperature from spectra recorded as a function of 
temperature or pressure are summarized. This section contains estimates 
of the dimensions of the mixing region, i.e., the spatial resolution of 
the measurements, and of the absolute magnitude of the scattered light. 

5.1.1 Reduction of Experimental Spectra for Comparison with Theory 

A typical single-shot CARS spectrum of room temperature N 2 as it 
was recorded by the OMA is shown in Fig. 5.1(a). The large, uneven back- 
ground is due to a combination of dark noise from the camera target and 
electronic noise associated with the readout and analog-to-digital con- 
version processes. The first correction applied to the raw data is the 
subtraction of this background to obtain a spectrum with a flat baseline. 
This is accomplished by recording a "no light" frame for which the OMA 
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scan parameters are identical to those used in recording the experimental 
spectrum; this no light frame is then subtracted from the signal frame. 

It is impractical to record a background frame for each signal frame, but 
because of constant drift in the background, it is necessary to record a 
new background correction frame every few minutes. The background- 
corrected spectrum corresponding to Fig. 5.1(a) is shown in Fig. 5.1(b). 
Although this correction appears to the eye to work extremely well, it 
does introduce baseline errors that lead directly to errors in the tempera- 
tures deduced in the fitting routine. For example, a single series of CARS 
spectra analyzed using two consecutively recorded background frames can 
yield average temperatures differing by more than the rms error. 

The second correction that must be applied to the data is one that 
compensates for the nonuniform response of the system as a function of 
position in the image plane. The major contributions to this response 
error are nonuniform throughput of the spectrometer and imaging optics 
and the nonuniform response of the detector itself. These errors are 
eliminated by dividing the background-corrected CARS spectra by a response 
frame consisting of a background-corrected white light spectrum obtained 
by placing the filament of a tungsten lamp at the sample region focus in 
the CARS set-up (see Fig. 3.3). White light from the lamp is directed 
along the same optical paths traversed by the sample and reference CARS 
signals, is dispersed by the spectrograph, and is imaged onto the OMA. 

The white light response frame therefore reflects nonuniformities in both 
the throughput of the entire optical system and the OMA response. 

Because the camera response is not strongly time or temperature dependent, 
a single response frame can be used with all data taken with a given 
optical alignment. A typical response spectrum is shown in Fig. 5.1(c). 
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Fig. 5.1. (a) Single-shot, uncorrected CARS spectrum as recorded 

by the OMA. (b) Same spectrum as shown in (a) , with 
background subtracted. (c) Typical white light spectrum 
used to correct the CARS spectra for nonuniform system 
response. 
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The integrated CARS power in a particular rotational line of a 
background- and response-corrected spectrum is obtained by summing 
the signals from all OMA channels in the vicinity of the line. The 
range of the summation should ideally be equal to the spacing between 
adjacent lines and should be centered on the line of interest. The 
full range of integration was employed for averaged spectra in which 
the baseline noise is small, but the range was reduced for single-shot 
spectra in order to avoid errors introduced by the noisy baseline 
between peaks where the signal is indistinguishable from the background. 
The final integrated-intensity distribution corresponding to a given 
experimental spectrum is obtained by dividing the integral of each CARS 
line by the total CARS signal from all lines. The resulting normalized 
distribution indicates what fraction of the total CARS intensity appears 
in each rotational line and can easily be compared to the calculated 
distributions . 

Spectra were recorded over a range of temperatures (135 to 296 K) 
at atmospheric pressure and over a range of pressures (100 Torr to 15 atm) 
at room temperature. Both single-shot and averaged spectra were recorded 
as a function of temperature, while all pressure-dependent spectra 
reported here are averaged over 10 or 15 laser firings. A series of 
5 or 10 spectra were recorded at each condition of temperature and 
pressure, and an average temperature and rms deviation were calculated 
for each series. Each experimental spectrum was compared with a series 
of theoretical distributions calculated at 1 K increments over the 
appropriate temperature range and at the appropriate pressure. The 
comparison was performed by calculating for each theoretical 
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distribution a sum of squares error parameter defined by 


e(T) 


l 

J 



cal 

J 


where is defined in Eq. (4.1). For room temperature spectra, the 
sum was usually over J = 4 to J = 14; at colder temperatures the range 
of the sum was reduced. The temperature of the experimental spectrum 
was then taken to be that for which e(T) was minimized. 

In order to fit all of the data, theoretical distributions were 
calculated from 100 to 330 K at atmospheric pressure, and from 270 to 
340 K at each 100 Torr increment in pressure over the range of 100 

I 

to 600 Torr and at higher pressures of 2.79, 4.57, 8.14, 11.71, and 
15.29 atm. Calculation of a complete set of distributions using the 
single resonance theory [Eq. (4.8)] required tens of hours of processor 
time on a mainframe computer system — hundreds of hours would have been 
required had the double resonance theory [Eq. (4.6)] been used. Once 
the distributions were available, however, an experimental spectrum 
could be background- and response-corrected and fit to yield a tempera- 
ture in several seconds using a microcomputer. 

5.1.2 Temperature and Intensity Measurements in N „ at STP 

A large number of spectra were recorded with the sample and 
reference gases at room temperature in order to determine the accuracy 
of the overall fitting procedure, including the background and response 
corrections, the reference normalization and the comparison with theory. 
Figure 5.2 shows the room temperature sample (b) and reference (a) spectra 
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Fig. 5.2. Single-shot spectra of room temperature sample gas (b) 

and reference gas (a). The distribution of signal among 
the rotational lines used to fit the temperature is shown 
in (c) for the reference-normalized sample spectrum and 
for the best-fit calculated spectrum. T ^ = 296 K, Tp, — 
293 K. 




obtained from a single laser firing. Neither resembles the calculated 
spectrum at 296 K [Fig. 4.2(a)] because of the spectral structure in the 
dye laser. Figure 5.2(c) shows the integrated intensity distributions 
corresponding to the reference-normalized sample spectrum and the best- 
fit calculated distribution. In this case, the thermocouple temperature 

was T = 296 K, the best-fit temperature was T_ = 293 K, and the average 

xC. r 

temperature of the series was <T > = 291.2 ± 4.0 K. Figure 5.3 is 

r 

similar to Fig. 5.2, only in this case the experimental spectra are 
averages of the spectra from 10 consecutive laser shots; again T^, = 296 
while T = 292 K and <T > = 291.1 ± 0.7 K. 

r r 

The distribution of temperatures obtained from 100 consecutive 
single-shot spectra is shown in Fig. 5.4(a); the average temperature 
for this series was 297.0 K and the nns error was ± 6.1 K, as indicated. 
The scatter in temperature may be due either to actual differences in 
the sample and reference spectra or due to noise introduced in 
recovering the signals from the OMA. In order to determine the relative 
importance of these two sources of error, additional spectra were 
recorded with the pump and dye beams blocked after the sample region. 

The sample CARS signal was then split, and half was directed along its 
usual path to the sample track of the OMA, while the other half 
followed the usual path of the reference signal to the reference track 
of the OMA. The two identical spectra recorded in this way were then 
processed as independent sample and reference spectra, and a best-fit 
temperature was obtained for each "sample" spectrum. The resulting 
temperature distribution for 50 consecutive shots is shown in Fig. 5.4(b) 
the series average was 295.6 ± 6.0 K. Since the sample and reference 
spectra were identical, the scatter can be attributed entirely to 
noise introduced by the OMA. The similar breadth of the distributions 
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in Figs. 5.4(a) and (b) suggests that this is the case even when 
separate sample and reference spectra are generated. 

In addition to correcting spectral errors in the sample gas 
spectrum, the reference normalization should also reduce the shot-to- 
shot intensity fluctuations in the sample CARS signals. The effective- 
ness of this amplitude normalization can be inferred from Fig. 5.5, 
which shows the distribution of the integrated intensity of 100 single- 
shot sample spectra in (a), and the ratio of sample to reference integrated 
intensities for the same 100 shots in (b) . For the unnormalized signals, 
the rms deviation from the mean intensity was ± 21%, while the rms 
deviation from mean of the normalized intensities was ± 5.8%. These 
results should be compared with those presented in Fig. 5.6, which shows 
the 50-shot distributions of the total sample signal intensity (a) and 
sample/reference intensity (b) obtained with a single spectrum divided 
into two parts, as above. The deviation from mean of the unnormalized 
signals was ± 24%, while the rms deviation of the normalized signals 
was ± 2.6%. Noise introduced by the OMA is responsible for the width 
of the distribution shown in Fig. 5.6(b), but the larger rms error 
observed in Fig. 5.5(b) indicates that there is significant inten- 
sity variation between the separate sample and reference signals 
apart from that introduced by the OMA. Thus, even with an ideal 
detector, attempts to determine the gas density from the normalized 
signal intensity will be limited to an accuracy of several percent; 
with the instrumentation used in this work, relative density measure- 
ments with better than ± 10% accuracy should be possible. 
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Fig. 5.3. Averaged spectra of room temperature sample gas (b) and 
reference gas (a) . The distribution of signal among the 
rotational lines used to fit the temperature is shown in 
(c) for the reference-normalized sample spectrum and for 
the best fit calculated spectrum. Tjc “ 296 K, Tp « 292 K; 
average is over 10 laser pulses. 
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Fig. 5.4. (a) Distribution of temperatures obtained from 100 single- 

shot sample and reference spectra; T = 296 K, <Tp> = 297.0 ± 
6.1 K. (b) Distribution of temperature obtained from 50 
consecutive single-shot sample spectra split into two parts; 
T = 296 K, <Tp> = 295.6 ± 6.0 K. Arrows in (a) and (b) 
indicate <Tp>, horizontal bars indicate rms deviations from 
the mean. 
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Fig. 5.5. Distributions of unnormalized (a) and reference-normalized 
(b) intensities of 100 consecutive single-shot CARS spectra 
of room temperature N 2 - rms deviations from the mean are 
± 21% for the unnormalized intensities and ± 5.8% for the 
normalized intensities. 
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Fig. 5.6. Distributions of unnormalized (a) and reference-normalized 
(b) intensities of 50 consecutive single-shot CARS spectra 
of room temperature N 2 . Sample and reference spectra were 
obtained by splitting the sample spectrum and are identical, 
rms deviations from the mean are ± 24% for the unnormalized 
intensities and ± 2.6% for the normalized intensities. 
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5.1.3 Temperature Measurements in Cold 

CARS spectra were recorded from ^ at atmospheric pressure over 

the temperature range of 135 to 296 K. Two series of ten spectra each 

were recorded at each temperature point; one series consisted of 

consecutive single-shot spectra, while the spectra in the other series 

were averaged over ten shots. Figure 5.7 shows the cold sample gas 

and the room temperature reference gas spectra obtained in a single 

shot. The integrated-peak intensity distribution for the reference- 

normalized sample gas spectrum and the best-fit calculated distribution 

are shown in Fig. 5.7(c). Only those peaks with good signal/noise 

in both the room temperature and cold spectra were used in the fitting 

procedure. For the spectra shown in Fig. 5.7, the thermocouple 

temperature was = 169 K, the fit temperature was T^ = 173 K, and 

the average temperature of the series was <T > = 170.2 + 5.4 K. 

r 

Figure 5.8 is similar to Fig. 5.7, except the sample and reference 
spectra were averaged over 10 shots; in this case T = 146 K, 

X *w< 

T = 140 K, and <T > = 141.8 ± 1.6 K. 

r r 

The results of the measurements in cold ^ gas are summarized in 
Fig. 5.9, which shows the CARS temperature plotted against the 
thermocouple temperature over the range 135 - 296 K. The error bars 
correspond to the rms deviation from mean of the single-shot measure- 
ments. The rms deviations from mean for the series of averaged 
spectra are not shown but are always smaller than the single-shot 
deviations. 

Figure 5.9 indicates very good agreement between the CARS and 
thermocouple temperatures, but it should be compared with Fig. 5.10, 
which shows the temperatures obtained from the same data used 
in Fig. 5.9, also plotted as a function of thermocouple temperature. 
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Fig. 5.7. Single-shot spectra of cold sample gas (b) and room 

temperature reference gas (a). The distribution of signal 
among the rotational lines used to fit the temperature 
is shown in (c) for the reference-normalized sample spectrum 
and for the best-fit calculated spectrum. Tjc = 169 K, 

Tp = 173 K. 
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Fig. 5.8. Averaged spectra of cold sample gas (b) and room temperature 
reference gas (a) . The distribution of signal among the 
rotational lines used to fit the temperature is shown in 
(c) for the reference-normalized sample spectrum and for 
the best-fit calculated spectrum. T TC = 146 K, T F = 140 K; 
average is over 10 laser pulses. 
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Comparison of CARS temperature with thermocouple temperature 
over the range 135 to 296 K. The average and nns deviation 
from mean are indicated for a series of 10 single-shot 
spectra at each temperature. The average temperature for 
a series of 10 averaged spectra is also shown, where each 
spectrum is averaged over 10 laser pulses. 
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5.10. Comparison of CARS and thermocouple temperatures for the 
same data as used in Fig. 6.9, but with an additional 
baseline correction used in the analysis. 



The only difference between the two figures is in the way the back- 
ground was subtracted from the raw spectra (see Fig. 5.1). The 
results in Fig. 5.9 were obtained by setting to zero the minimum 
of each spectrum once the background had been subtracted; the signal 
from each line was then integrated according to the procedure out- 
lined earlier in this chapter. The correction used to obtain the 
results of Fig. 5.10 was similar, but in addition to setting the 
minimum of each background-subtracted spectrum to zero, an individual 
baseline was also subtracted from each peak, essentially setting 
to zero the local minimum in the vicinity of each rotational line. 
Although this step seems somewhat artificial, experience has shown 
that in many cases it does help to reduce the scatter in the temperature 
fits, and occasionally improves the accuracy at room temperature. 
Unfortunately, neither correction is successful in all cases, so it 
appears that the irreproducibility of the background of successive 
OMA frames is the principal limit to the accuracy of this temperature 
measurement technique. 

The most important source of error in these measurements is the 
OMA. The disagreement between the thermocouple and CARS temperatures 
shown in Fig. 5.10 appears to be due to baseline errors caused by 
the inexact background subtraction on the OMA rather than to real 
spectroscopic effects or to inadequacies of the theory used for the 
calculation of the theoretical distributions. The scatter in the 
temperature measurements can also be attributed to noise introduced 
by the OMA; with a better detector there would probably still be some 
scatter due to real differences in the sample and reference spectra, 
but it is not possible to conclude from the present work what the 


ultimate limit would be. 
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5.1.4 Temperature as a Function of Np Pressure 

Room Temperature spectra were recorded over a wide range of 
pressures in order to illustrate the importance of the interferences 
between the individual resonant Raman lines. At low pressure, the 
pressure-broadened lines are very narrow (^.005 cm - ^ at 100 torr) 
and the peak resonant susceptibility is very high ('Vf ), so the 
effect of the nonresonant background due to and the tails of other 
off-resonant lines is small. In this case the intensity distributions 
are very similar to those obtained assuming isolated lines with no 
X NR . As the pressure, and therefore the linewidths, increase, 
however, the peak resonant susceptibility decreases and the influence 
of the real background becomes more pronounced. The effect on the 
integrated- intensity distributions is an increase in magnitude of the 
low J lines, since they have the broadest linewidths and are thus most 
susceptible to interference effects. The high J lines are also 
enhanced because they interact with the off-resonant (negative value ) 
tails due to the lower level resonances, just as the low J lines see 
the sum of all of the off-resonant (positive value) tails due to 
the higher J resonance^. The enhancement of the high J lines is 
less than that of the low J lines, however, because of the much 
narrower linewidths of the higher rotational resonances. The over- 
all effect of these interferences is a shift toward lower J levels 
in the integrated CARS intensity distribution with increasing 
pressure at constant temperature. Thus, if spectra from two samples 
at the same temperature but different pressure are compared using 
the same calculated distributions, the higher pressure spectrum 
will appear to be colder than the lower pressure spectrum. 
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The experimental results that demonstrate this apparent change 
of temperature with pressure are summarized in Fig. 5.11, where the 
CARS temperature is plotted as a function of pressure. Each point 
is the average temperature of a series of 10 spectra, each of which 
is averaged over ten consecutive laser pulses. The spectra were 
recorded with the reference gas at 1 atm and 296 K, while the sample 
gas pressure was varied at constant temperature (296 K) . The solid 
circles and squares represent the temperatures obtained using 
theoretical intensity distributions calculated for the appropriate 
sample gas pressure, while the open circles and squares represent 
the temperatures obtained by analyzing all the data with distributions 
calculated for 1 atm pressure. The 6 K offset at 1 atm is an 
experimental artifact — the low and high pressure spectra were recorded 
at different times, and the response correction used with the high 
pressure results systematically weighted the sample gas spectra 
toward higher temperature. The conclusion to be drawn from Fig. 5.11 
is that is is necessary to know the Raman linewidths of the ^ sample 
gas in order to accurately determine the temperature from its CARS 
spectrum. 

The measurements of the temperature as a function of pressure are 
a good test of the effectiveness of the resonant reference normalization. 
The fundamental difficulty with the use of a resonant reference is that 
the normalization is not exact if the sample and reference gases have 
different Raman linewidths. Although the integration in Eq. (4.2) was 
performed by assuming that the dye laser spectrum is flat, it was 
recognized that the CARS process itself performs a similar integration 
over the exact structures that do exist in the dye laser spectrum. As 
long as the sample and reference gases have the same linewidths, the 
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Fig. 5.11. CARS temperature as a function of sample gas pressure for 
room temperature sample and STP reference. Solid squares 
and dots indicate average temperatures obtained using 
theoretical distributions calculated for sample gas 
pressure; open squares and dots indicate average temperatures 
obtained using theoretical distributions calculated for 
1 atm pressure. Error bars show the rms deviation from 
mean for a series of 10 averaged spectra at each pressure. 



results of this physical integration will be the same in both media. 

If the linewidths are different, however, the integrand ( [ x ^ \ ) in 

Eq. (4.2) weights the same dye laser spectral component differently 
in the sample and reference spectra, leading to slightly different 
values for the two integrated signal strengths. The error introduced 
in this way should become more important as the difference between 
the pressures of the sample and reference gases increases. For example, 
the signal from a very broad rotational line may be significantly 
influenced by dye laser spectral components that are too far away to 
be of importance to the signal from a narrower line. However, this 
sort of error is not observed in this work; in fact, the results 
shown in Fig. 5.11 indicate that the normalization is successful even 
when the linewidths differ by a factor of 15. This may be attributed 
to the broad linewidth of the pump laser used in these experiments 
since even at a pressure of 15 atm , the Raman linewidths are still 
less than the laser linewidth. The CARS linewidths are therefore 
determined largely by the pump laser linewidth and, as a result, are 
relatively independent of pressure. The use of a pump laser with 
a large linewidth is therefore an advantage if measurements are to be 
made over a wide range of pressures. 

Two new problems are encountered at pressures below 100 torr. The 
first is that the contribution to the Raman linewidths due to Doppler 
broadening becomes important, thus altering the shapes and pressure 
dependence of the rotational linewidths. The analysis presented in this 
work assumed Lorentzian lineshapes appropriate to homogeneous broaden- 
ing and cannot be easily patched up to deal with the Voigt 
lineshapes appropriate when the Doppler and pressure broadened 
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linewidths are comparable. For the experimental geometry used in this 
work, the Doppler linewidth at room temperature is 0012 cm * ; the 
pressure-broadened linewidth of the J=10 line reaches this value 
at ^10 Torr so 100 Torr represents a safe lower limit to the range 
of pressures in which Doppler broadening can be neglected. 

The second problem at low pressure is that the assumptions 
implicit in the theory used to calculate the rotational spectra break 
down. In evaluating the integrals in Eq. (4.2), it was assumed that 
the laser intensity is distributed continuously with frequency over 
the laser linewidths. This is not the case, however, since the laser 
outputs consist of many modes spaced ^.003 cm ^ apart in the YAG laser 
output and ^.025 cm * apart in the dye laser output. Equation (4.2) 
should be valid only in the limit that the Raman linewidths are broad 
compared to the laser mode spacing; this is true only for the YAG laser 
at atmospheric pressure (T^ ^.05 cm *) and is true for neither laser 
at 100 torr (r ''-.006 cm *). In fact, the results shown in Fig. 5.11 
do indicate that the temperature fitting routine is more successful 
with the sample pressure increased by a factor of 15 than it is with 
the pressure decreased by only a factor of 7.5. 

5.1.5 Spatial Resolution 

The spatial resolution of a CARS measurement is determined by 
the dimensions of the region in which the nonlinear mixing process 
takes place. In experiments employing colinear phase-matching, this 
interaction region is roughly cylindrical in shape; its length is 
determined principally by the speed of the focusing lens since the 
mixing efficiency is greatest in the region of tightest focusing. 

In experiments that use noncolinear phase-matching (BOXCARS, 
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three-dimensional phase-matching), the interaction region is also roughly 
cylindrical, but the length is determined by the phase-matching geometry 
since a signal can only be generated from the region of overlap of the 
three input laser beams. The transverse dimensions of the scattering 
region are difficult to measure directly, but they are generally deter- 
mined by the focal waists of the laser beams which are on the order of 
100 pm in diameter. The longitudinal dimension [%£ in Eq. (1.8)] can 
be measured directly and is usually much greater than the transverse 
dimensions, so that, in practice, the actual resolution is limited by 
the interaction length. 

The interaction length was determined for the geometry used in 
this work by measuring the nonresonant CARS signal generated in a 
thin piece of microscope cover glass as a function of the longitudinal 
position of the glass in the focal region. The OMA was replaced by 
a photomultiplier tube in making this measurement because the 
nonresonant signal from the glass was quite small (the laser inten- 
sities had to be greatly reduced to avoid damage to the glass). 

The detected CARS intensity as a function of position z generated 
in a glass slide of thickness Az is approximately 

l| lasS (z) - ] X | 2 |(eJe*) z | 2 Az 2 , ( 5 . 1 ) 

2 * 

where (EE) is the product of the input laser fields at position z. 

L b Z 

In the gaseous experiment the signal is generated from the entire 
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region and is given by 


r gas 


X 0&), dz) 2 


% I I X (e£e*) z Az ! Z ’ 


(5.2) 


where the integration and summation are over the entire interaction 
region. Equations (5.1) and (5.2) can be combined to yield 

lf s - [ I (lf aSS (z)) 1/2 ] 2 . 

z 

The best estimate of the interaction length, therefore, is given not 

slsss 

by the distribution of if with longitudinal position z, but by 

. 4 1/2 
R J.3S S 

the distribution of (I| ) with z. Failure to recognize this 

point results in an underestimate of the interaction length. The 
above analysis yields an interaction length of 2.5 to 3.0 mm for 
the geometry used in this work. 

5.1.6 CARS Signal Magnitude 

In order to make an estimate of the magnitude of the scattered 
CARS power, it was necessary to make an absolute calibration of the 
response of the OMA to incident light. This was accomplished by 
using a photodiode of known responsivity at 532 nm to determine the 
average energy in a series of greatly attenuated pulses from the YAG 
laser. This number for the average energy was then used in conjunc- 
tion with the integrated area of the laser line detected with the 
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OMA to determine a calibration factor; for the scan format used 
in most of this work, the calibration factor was 'v. 22 photons /count , 
where 1 count corresponds to the least significant bit of the OMA's 
analog-to-digital converter . 

Spectra of room temperature and pressure could be recorded 
with good signal/noise with a total laser energy at 532 nm of 20 mj/ 
pulse and a total broad-band dye laser energy of 1.0 mJ/pulse. 
Application of the calibration factor of 22 photons/count to the 
strongest rotational line (J = 10 ) in a typical spectrum yields a 
total signal at the detector of "v 10 ^ photons for the single line. 

The signal expected from an individual CARS line can be easily 
calculated using Eqs. (1.18) and (1.8). The following values are 
appropriate to the above experimental conditions: 


( 37 - ) ^ 10( 77 ) = 3.8 x 10 30 cm 2 /mol - sr 

dQ rot dn Q 


A = .033 (J = 10) 

T = .038 cm ^ 

19 3 

N = 2.4 x 10 mol/cm 
L = 1.5 mm 

2 

I = 10 mJ /10 nsec/r (100 pm) 

L 

2 -1 

1=1 mJ/10 nsec/r (100 pm) /150 cm bandwidth. 

l) 


Because the signal from only a single line is desired, the integra- 
tion over cd can be performed analytically to yield 
A 
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9 


if interferences with x NR and other rotational lines are neglected. 

This calculation yields a scattered signal from the J = 10 line of 
9 

'v* 10 photons. This must be multiplied by a factor of ^ .15 to 
account for the throughput of the optics and spectrograph in order 

g 

to obtain an estimate of ^ 2 x 10 photons at the detector. This 
estimate is 200 times greater than the observed signal, but dis- 
crepancies on this order are common in pulsed CARS experiments 
due to the nonlinear nature of the mixing process and to the non- 
ideal temporal and spatial properties of the laser beams. 

5.1.7 Discussion 

The experimental results presented in Section 5.1 indicate that the 
broad-band, single-pulse CARS technique is capable of accurate measure- 
ment of the rotational temperature of ^ gas at room temperature and 
below. Because of the broad-band approach, discussion is made of the 
novel usage of a resonant reference cell to normalize the wavelength 
dependent dye laser intensity and of the problems associated with the 
OMA. The good agreement between the thermocouple and CARS temperatures 
shown in Fig. 5.9 suggests that the theory used in calculating the 
integrated intensity distributions is adequate for this application 
and indicates that the change in J-dependence of the rotational line- 
widths with temperature does not represent a significant source of 
error. It should be possible, therefore, to extend this technique to 
the measurement of rotational temperature in other molecules, since the 
theory used in analyzing the data in the current work is not unique 
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to N^. It should also be a straightforward matter to extend the 
temperature range over which this technique is applied to both colder 
and warmer temperatures, although the analysis will have to be 
slightly modified in the latter case to include the effects of the 
rotational-vibrational interaction and the contribution of the 
rotational hot bands. 

A natural extension of this work would be to the simultaneous 

measurement of pressure as well as temperature in the scattering 

medium. Just as the temperature of the sample gas is deduced by 

comparing the spectral features of the sample and reference spectra, 

the sample gas pressure can be deduced by comparing the intensities 

of the two spectra. Unfortunately, the sample gas pressure cannot be 

obtained directly form the ratio of the signal intensities, because 

the overall amplitude of a rotational CARS spectrum depends strongly 

on temperature as well as on pressure. The temperature dependence 

(3) 

comes about through three terms in x : the molecular density, which 
increases with decreasing temperature (at constant pressure); the 
Raman linewidths, which increase with decreasing temperature; and the 
population difference factor, which increases with decreasing temp- 
erature due to the condensation of the rotational population into the 
low J levels. The overall effect is an increase in the magnitude of 
the CARS signal with decreasing temperature; for example, the integral 
of the strongest rotational line at 150 K (J=6) is 10.7 times greater 
than that of the strongest line at 300 K (J=10). This presents a 
problem for the determination of temperature and pressure from a single 
pulse because it is necessary to know the gas temperature in order 
to obtain the pressure, but we have seen that accurate temperature 
fits require knowledge of the pressure. One possible solution to this 
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dilemma is to solve for the temperature and pressure by an iterative 
technique in which an initial estimate is made for the temperature, 
assuming a pressure of 1 atm. A pressure based on this temperature can 
then be deduced and used to refine the temperature fit, which will 
allow further refinement of the pressure estimate. An analysis of this 
type could be applied to the spectra recorded in this work to yield 
time-resolved measurements of pressure and temperature with good spatial 
resolution and perhaps even to multipoint determination. 

5.2 Rotational and Vibrational CARS Measurements of in a Flame 

Coherent anti-Stokes Raman scattering (CARS) has been used exten- 
sively as a diagnostic technique for temperature measurement in combus- 
tion systems (Regnier and Taran, 1973; Stenhouse et al., 1979; Rahn 
et al., 1979; Hall and Eckbreth, 1981; Eckbreth and Schreiber, 1981). 

The relative intensities of the vibrational-rotational bands of Nj gas 
provide a sensitive temperature probe at flame temperatures. Pure 
rotational CARS is superior to vibrational CARS for thermometry of ^ 
gas at room temperature and below, since the spectra are easily 
resolvable in contrast to the congestion of the rotational lines in the 
vibrational bands of the Q-branch spectra. Furthermore, the peak third- 
order susceptibility for rotational CARS is larger than that for vibra- 
tional CARS (Fenner et al., 1973; Penney et al., 1974; Kondilenko et al., 
1980). However, at higher temperatures, rotational CARS signal intensi- 
ties are reduced substantially as the population difference factor 
between the initial and final states for the CARS process approaches 
zero. In order for this population difference to become significant, 
the initial states need to be of high J, and thus are less populated. 

This small population difference at higher temperatures has led most 
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workers to assume that pure rotational CARS has much smaller signal 
strengths relative to vibrational CARS for systems at flame tempera- 
tures. We investigated this assumption by comparing, under similar 
conditions in the same laboratory, the rotational and vibrational CARS 
techniques for the determination of the Nj temperature in a flame around 
1600-1900 K. 

The flame was generated with a flat-flame burner and a premixed 
flow of CH^ and air was used. The thermocouple was fabricated from a 
0.05 mm diameter Pt and Pt~13%Rh wire yielding a junction bead diameter 
of 0.15 mm (Farrow et al., 1982). As with the room and low temperature 
CARS experiments the second harmonic (532 nm) of a Nd:YAG laser was 
used as the monochromatic pump radiation. For the rotational CARS, 
the third harmonic of Nd:YAG (355 nm) was used to pump a coumarin 540A 
broadband dye laser (532-542 nm) . The frequency selection element for 
this dye laser configuration was an angle-tuned dielectric mirror 
coated for total reflection at 532 nm, 45° incidence and s-polarization. 
This element provided a very smooth tuning curve and allowed lasing at 
the red edge of the gain curve, a characteristic the original element 
did not possess. For the vibrational CARS, a portion of the second 
harmonic of Nd:YAG (532 nm) was used to pump a broadband rhodamine 610 
broadband dye laser (600-610 nm) . The dye laser consisted of a trans- 
versely pumped oscillator and amplifier, with frequency selection 
achieved by angle-tuning the original dielectric filter placed in the 
oscillator cavity. 

The small-angle three-dimensional phase-matching geometry was 
used for both the rotational and vibrational CARS. The spectrograph 
used for this work, as well as for the remainder of the experiments 
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presented in this report, was a 0.5 m spectrograph with vertical slits 
and a 2400 line /mm holographically ruled grating. A reference cell of 
Ar at 9.5 atm was used to generate a nonresonant CARS signal which 
essentially mimicked the spectral intensity fluctuations of the coumarin 
and rhodamine dye laser, and which occured in the rotational and 
vibrational-rotational CARS spectrum, respectively. This reference 
signal was imaged on the lower portion of the spectrometer slits. The 
CARS of N£ and the reference were detected simultaneously on spatially 
distinct portions of the SIT vidicon. The overall spectral resolution 
of the spectrograph and vidicon system was 2.7 cm The flame and 
reference signals were corrected individually for the vidicon background 
and for the overall system response. The corrected flame signal was 
normalized by the corrected nonresonant Ar reference signal. 

A typical vibrational CARS spectrum is presented in Fig. 5.12, with 
the experimental and the calculated spectra represented by the solid and 
dashed lines, respectively. For the vibrational CARS spectra, the 
calculated best least-squares fit temperature of 7 spectra was 1841 K 
and the radiation-corrected thermocouple measurement was 1865 K. The 
uncertainty in the vibrational CARS deduced temperature is ^3%. The 
vibrational temperatures were determined by Dr. Roger L. Farrow of 
Sandia National Laboratories, Livermore, California. These temperatures 
were deduced from nonlinear least squares fits of theoretical to experi- 
mental spectra. The theoretical spectra were obtained by calculating 
the complex Raman susceptibility, adding the appropriate nonresonant 
susceptibility and convoluting the resulting squared modulus with an 
instrumental lineshape function. These calculations are described more 
fully elsewhere (Farrow et al., 1982). 
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Fig. 5.12. Vibrational CARS spectra of the N 2 Q-branch, where the 
solid line is the experimental spectrum and the dashed 
line is the theoretical spectrum. The thermocouple 
temperature was 1865 K and the best least-squares fit 
CARS temperature was 1864 K. 
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A typical pure rotational CARS spectrum of Ng i n a flame is presented 
in Fig. 5.13. The experimental spectrum is shown in Fig. 5.13(a) and a 
calculated spectrum, chosen for best visual fit, is shown in Fig. 5.13(b). 
The maximum intensity is for around J = 24. The method of determining 
the temperature from the spectral intensity distribution is different 
for the two CARS techniques. The approach taken to deduce the rotational 
temperature is to calculate the relative integrated peak area of all the 
J peaks over a sufficiently wide range of temperature and then to deter- 
mine the least squares difference between the experimental and the 
theoretical integrated intensity of all the J peaks. This analysis would 
be similar to that described for the rotational CARS of room temperature 
and cold Nj gas except that the influence of hot bands should be taken 
into account in calculating the intensity distribution of the J peaks. 

We are currently working on the calculation of the theoretical intensity 
distribution for temperatures up to 2000 K. 

The higher efficiency of the second-harmonic/rhodamine dye laser 
system relative to the third-harmonic/coumarin dye laser system has 
been cited as a primary reason for choosing vibrational CARS over rota- 
tional CARS (Hall and Eckbreth, 1981). For the vibrational CARS, the 

+ was 120 mJ and the energy of the broadband 
rhodamine dye laser <uig> was 3.0 mJ /pulse. For the rotational CARS, 
the total pump energy was only 100 mJ, since portions of the fundamental 
and second harmonic were converted to the third harmonic in order to 
pump the coumarin dye laser . The broadband coumarin dye laser energy 
was 2.6 mJ /pulse. Since the magnitude of the CARS signal is proportional 
to the triple intensity product, I (co^) I (ti£) I (ojg) , for our experimental 


total pump energy (oj^ 
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Fig. 5.13. Typical pure rotational spectrum of ^ gas in a flame. 

(a) Experimental spectrum and (b) calculated spectrum 
chosen for best visual fit. 
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situation, the vibrational CARS signal had an input intensity advantage 

of n.l.5x, all other factors being equal. 

The relative signal strengths for vibrational and rotational CARS 

have been calculated by Hall and Eckbreth (1981) . Assuming the same 

triple intensity product, with the spectral linewidth equal to 0.8 cm * 

and and equal to 150 cm * for <a)g>, the ratio of vibrational 

CARS intensity to rotational CARS intensity was calculated. At T = 300 K, 

I ., / 1 . * 1.25 and at T = 1800 K, I / 1 . = 2.5. Our measurement 

vib rot vib rot 

of relative signal strengths did not support this estimation. We 

measured signal strength for both the rotational and vibrational CARS 

by introducing calibrated neutral density filters into the CARS beam 

until their signals were indistinguishable from the vidicon readout noise. 

The measured rotational and vibrational CARS peak signal strengths under 

the same flame conditions CVL800 K) were nearly the same, although the 

rotational CARS seemed somewhat stronger. Since the triple intensity 

product for vibrational CARS was nd..5X of that for the rotational CARS, 

we conclude that the intensity ratio should be I / 1 ^ < 1. That is, 

vib rot 

for the same intensity product, the peak I is 1-2X more intense the 
1 even at 1800 K. 

The sensitivity of the rotational and vibrational CARS spectra to 
temperature changes must be considered in addition to total signal 
intensity. Rotational CARS has the advantage of not requiring higher 
spectral resolution than ^3 cm ^ since the peaks are adequately separated 
and only the integrated peak area for each J is required. The vibrational 
CARS procedure requires the fitting of the lineshape and higher resolution 
does improve the accuracy of the temperature deductions (Klick et al., 1981). 
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The main sources of inaccuracy encountered with rotational CARS 
measurements are associated with the broadband dye laser. The first 
is the decreased output of the coumarin dye laser relative to that of 
the rhodamine laser, and particularly the necessity to pump the coumarin 
dye with the third harmonic, which decreases the second-harmonic intensity 
which is used for the two pump beams. The second disadvantage is that any 
spectral structure in the dye laser output tends to influence the intensi- 
ties of the isolated, relatively narrow rotational lines (FWHM ^2 cm *) 
much more critically than those of the broader vibrational band. Because 
of distortions in the optical system (spectrograph and vidicon) it is 
not possible to obtain an accurate correspondence between the signal and 
reference channels with respect to spectral frequency. This difficulty 
was minimized by imaging the sample and reference signal on adjacent 
portions of the vidicon for which the spatial distortion was minimal. 

Even so, for the rotational spectrum, the peaks are on the order of a 
few channels wide, so even a single channel mismatch creates a significant 
error in the integrated area for various J modes. 

Nevertheless, in spite of several disadvantages associated with 
rotational CARS, it is still a viable technique for the measurement of 
flame temperatures. If the difficulties with the normalization could be 
solved adequately, rotational CARS would be particularly advantageous 
for systems requiring the measurement of a range of temperatures 
(77 K-1900 K) with the same experimental apparatus. 

5.3 Simultaneous Observations of CSRS and CARS 

Coherent anti-Stokes Raman spectroscopy (CARS) and coherent Stokes 
Raman spectroscopy (CSRS) are closely related nonlinear optical processes. 

3 

The form of the third-order susceptibility (X ) for CARS and CSRS differs 
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when the excited electronic states are nearly resonant with the photon 

3 

energy of either the pump wave or the probe wave. Depending on the x for 

the nonresonant process, different line shapes for the CSRS and CARS 

spectra can result. In fact, such line shape behavior can be used to probe 

the dipole matrix elements between the ground state and a particular 

3 

excited state and to elucidate the characteristics of the x near elec- 
tronic resonance conditions (Oudar and Shen, 1980; Carreira et al., 1978; 
Igarashi et al., 1980; Dutta and Spiro, 1978). 

When the energies of the excited electronic states are far removed 
from the photon energies of any of the input beams, the CSRS and CARS 
processes depend on molecular parameters in precisely the same way 
(Druet and Taran, 1981) . Consequently, the CSRS and CARS spectra contain 
identical Information about the rotational and/or vibrational-rotational 
levels, their linewidths and the population difference factors between 
the initial and final levels. The CARS technique is generally preferred 
because of its excellent rejection of unwanted fluorescence and because 
It is more convenient to produce a Stokes, rather than an anti-Stokes, 
probe beam since a red-shifted dye laser can be pumped with the same 
radiation that provides the pump beam. 

We simultaneously observed the CARS and CSRS of the rotational levels 
of air and of ^ by using broadband generation and detection. Previously, 
spectra of pure rotational CARS and/or CSRS have been obtained sequentially 
by scanning the frequency of the probe radiation. The probe radiation was 
provided by a broadband dye laser with a spectral output which spans the 
Stokes and anti-Stokes regions. The CSRS and CARS signals were dispersed 
by a spectrograph and detected simultaneously by the OMA. The small-angle 
three-dimensional phase-matching geometry was used (a * 1.5° and 8 = Y = 1.5° 
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for the small shifts and angles involved* these angles are essentially 
equal). This configuration provided spatial resolution at the inter- 
action region and spatial separation of the CSRS and CARS beams from the 
input beams. For large Raman shifts, the three beam angles must be 
adjusted so that the phase-matching condition is satisfied either for 
optimum Stokes production for the CSRS process or for optimum anti-Stokes 
production for the CARS process. However, for small Raman shifts such 
as those encountered in pure rotational CSRS and CARS from ^ and C^, the 
phase-matching condition can be nearly satisfied for both CSRS and CARS 
processes. 

Essentially the same experimental arrangement was used as in the 
previous CARS work. The coumarin 500 dye laser was tuned so that the 
broadband radiation (^200 cm *) was centered at co^, i.e., the dye laser 
output provided both the broadband probe Stokes radiation, < to g > , and the 
anti-Stokes radiation, <oo A >. This particular phase-matching configuration 
provided adequate spatial separation of the broadband signals generated 
at <aj CSRS > and <U3 cars > f rom the three input beams and '“W- The 

= kg are shown in 

Fig. 5.14. The spectrally dispersed CARS and CSRS signals were detected 
simultaneously by the OMA. A small rectangular mask was inserted at the 
exit focal plane of the spectrograph in order to physically block the 

. The stray radiation at <l V and at <o)g> was much 
weaker than either the stray pump radiation or the CSRS and CARS signals. 
Thus, the three-dimensional geometry employed and the insertion of the 
mask provided sufficient rejection of unwanted radiation such that the 
CSRS and CARS signals could be observed simultaneously with the same 
detector, even though the stray pump radiation was not totally suppressed. 


stray radiation at ca^ 


/ 


relative directions of these beams, k^, k^, and k^ 
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k * 


k AS + k CSRS 


Fig. 5.14. (a) Three-dimensional phase-matching configuration for 

CARS process and (b) corresponding phase-matching for 
the CSRS process. The phase -mismatch for the CSRS 
process is indicated by Ak. 




The simultaneous CSRS and CARS spectra of room temperature air and 
^ are shown in Figs. 5.15(a) and (b) , respectively. These spectra were 
corrected for the vidicon background but were not normalized with respect 
to the variation of the dye laser intensity with frequency. The CSRS and 
CARS spectra are essentially mirror images. The three input beams were 
aligned to optimize the CARS process [see' Fig. 5.14(a)]. The slight 
variations in the intensity ratios of the CSRS and CARS peaks are due 
to the spectral structures inherent in the broadband dye laser which do 
not appear to average to a smooth distribution even after many pulses 
are accumulated on the vidicon target. 

The magnitude of the phase mismatch Ak depends on the phase-matching 
angles a and 3 selected for the three input beams, as well as on the 
frequency of the Raman shift, oi^. Figure 5.16 illustrates the dependence 
of the CSRS coherence length (l co j 1 = ir/Ak) on the phase-matching angles 
when the CSRS signal is detected at y» the optimum angle for the CARS 
signal. It is convenient to compare the coherence length for the CSRS 
signal to the interaction length (1 ) defined as the region in which 

the three input beams overlap. If all three beams have a diameter, d, 
then 1 can be approximated by d/sin2a. As long as 1 cq ^ > l^ nt » then 
the CSRS signal should be nearly the same strength as the CARS signal, 
since the latter also is limited by 1 Note that for large 0 ^ (e.g., 

the vibrational frequencies of Nj, C^, or Hj) 1^^ rapidly approaches 
1^ even for small-angle phase-matching configurations. For large u^, 
lcoh"* w ^ en a an ^ £ are greater than 5°. Consequently, a small-angle 

geometry is required for the simultaneous generation of CSRS and CARS 
for vibrational shifts. In practice, it would be difficult to generate 
the broadband Stokes and anti-Stokes probe radiation simultaneously. 
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PURE NITROGEN 
N 2 ( J J(yj ± 2 ) 
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CSRS 



RAMAN SHIFT 
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Fig. 5.15. Simultaneously observed experimental rotational CSRS and 
CARS spectrum (a) of ^ gas and (b) of air. 
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COHERENCE LENGTH (cm) FOR CSRS 


Fig. 


10 



PHASE- MATCHING ANGLE 01 FOR CARS 


5.16. Coherence length for the CSRS process as a function of the 
phase-matching angle a for the CARS process for various 
Raman frequencies. The interaction length as a function 
of a is indicated by the dashed line. 



In addition, the simultaneous observation with a single detector is 
problematic because of the large spectral range spanned in the vibra- 
tional case. However, it is easy to generate and observe both CSRS and 
CARS for the rotational shifts in molecules like Nj and C >2 since is 
small, even for rotational shifts from J = 0 to 15. Although the above 
discussion assumes a dispersionless medium, similar conclusions can be 
reached for the simultaneous generation of CSRS and CARS for liquid and 
solid media possessing small oj^. 

The small-angle phase-matching configuration used in this experiment 
is sufficient to exclude the detection of signals generated from stimu- 
lated Raman gain (Owyoung, 1978) since the output beams are not 
traveling in the same direction as the input probe beams. The inadvertent 
detection of stimulated Raman scattering due to the intense pump beam is 
also excluded since the CSRS and CARS signals disappear when the weaker 
^S+A 5 * P ro ^ e beam is blocked. 

Simultaneous measurement of nonresonant CSRS and CARS does not 
provide new information regarding the molecular system. However, similar 
experiments under resonant conditions may provide a means of investi- 
gating the possible interaction of the CSRS and CARS processes in driving 
the rotational levels. 
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PHASE-MATCHING FOR MULTIPOINT CARS 


The determination of instantaneous multipoint scalar parameters 
(temperature and/or density) in a turbulent flow has been one of the 
goals of this research. In particular, we have investigated the feasi- 
bility of using CARS to obtain instantaneous temperature and concentra- 
tion measurements from many spatially resolved points along a line 
intersecting the flow. In an essentially dispersionless medium like a 
dilute gas, there are various configurations which satisfy the phase- 
matching conditions discussed in Sections 1.1.1 and 2. There are several 
criteria to consider in choosing the optimal geometry: 1) the energy 

density of the input beams in the interaction region; 2) the spatial 
resolution of each volume element ; 3) the spatial separation of the 
anti-Stokes beam from the input beams; and 4) the coherence length for 
the frequency shifts of the rotational lines involved. In order to 
obtain spectral information from points along a line, the input beams 
must intersect in an elongated cylindrical volume such that there are 
several well-defined volume elements. Our goal was to have 10-20 
resolvable volume elements with dimensions approximately 0.1 mm x 
0.1 mm x 0.1 mm. 

6.1 Comparison of Various Multipoint Geometries 

The phase-matching schemes considered are shown in Fig. 6.1. Since 
the intersection of the beams occurs along a line, one or more of the 
input beams must be focused into a sheet of radiation using a cylindrical 
lens. The configuration shown in Fig. 6.1(c) is the same as the small- 
angle geometry used for the single point experiments. This has the 
disadvantage of requiring that all three input beams be focused into 
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( a )p!anar 

large-angle 


( b)pianar 

counter- propagating 




( c )folded small-angle (d)folded large-angle 




Fig. 6.1. Various phase-matching geometries for spatially-resolved 
multipoint CARS. The incoming pump beams are represented 
by £, and £ L < The Stokes beam is represented by kg and 
the generated anti-Stokes beam is represented by k^,. 


107 







sheets of radiation. This substantially decreases the available energy 
density in the interaction region. In addition, this geometry has very 
poor spatial resolution. The configuration shown in Fig. 6.1(a) is a 
planar, two-dimensional geometry. The spatial resolution for this con- 
figuration is essentially optimal, since the sheet of pump radiation, k^, 
is nearly perpendicular to both and lc . This arrangement requires 
only two of the input beams, ^ and to be focused into a sheet. The 
spatial discrimination of from is adequate for the relatively 
large frequency shifts which occur in vibrational CARS (Murphy et al., 
1979). However, in rotational CARS, where the shifts are small, the 
direction of is almost the same as that of and adequate spatial 
discrimination is not possible. A counter-propagating configuration, 
shown in Fig. 6.1(b), has the advantage of only requiring one of the pump 
beams to be focused into a sheet of radiation, considerably increasing 
the available energy density in the interaction region. This advantage 
is offset by very poor spatial resolution. The large-angle nonplanar 
three-dimensional geometry shown in Fig. 6.1(d) is the arrangement that 
was chosen for this work. This geometry, where ^ is exactly perpen- 
dicular to ^ and leg, has the optimum spatial resolution. Initially, 
it was thought that two beams (k and k 0 ) would have to be formed into 
sheets. It was determined that only needed to be focused into a 
sheet and that ^ and could be focused to a line with spherical lenses. 
The angle 6 between ^ and was small enough that there was sufficient 
overlap of these beams in the cylindrical volume defined by their inter- 
section with the perpendicular sheet formed by 5^. The anticipated 
difficulties with this configuration were the intensity of the CARS signal 
due to the extremely small interaction region (about 30X less than in the 
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small-angle single-point geometry) and the problems of achieving the 
necessary overlap and proper alignment of the three input beams. 

6.2 Sensitivity of Phase-matching Angles 

For the large-angle three-dimensional phase-matching geometry chosen, 
the selection of the angle a as 90° uniquely determines the remainder of 
the angles for a given frequency shift. It follows from the phase- 
matching condition that 


cos 6 58 - 2 caJj. 


k L 


6 = 5.671 


sin y = k£sin6 


Y = 5.616 


sin e = k - k/cos6 

_L h 

Ic cosy 
S 


e = 0.558 


where ol = 91.5 cm ^ is the frequency shift of the maximum J number 
(J = 10) at room temperature. 

The precision of the optical alignment is extremely critical to the 
success of these noncolinear CARS experiments. This is particularly true 
in the large-angle case where the interaction region is small. It was 
felt that it would be useful to know which of the angular adjustments 
were the most sensitive with respect to the phase mismatch Ak and to the 
coherence length, 1^^ * Tr/Ak. These quantities were calculated for it 
for angular displacements in the yz plane and in the xy plane (see 
Fig. 6.2). An angular displacement of ±0.5° in the yz plane (i.e., nearly 
perpendicular to the direction of gives Ak = 16 cm * and 1 ^ = 0.20 cm 
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An angular displacement of ±0.5° In the xy plane (i.e., nearly parallel to 
the direction of jc g ) gives Ak * 160 cm ^ and 1 cqR * 0.02 cm. The latter 
adjustment is 10X more critical than the former. To optimize the overlap 
of the input beams, it is clear that, given a choice, it would be more 
advantageous to change k R in a direction perpendicular to It g rather than 
parallel to 

In order to determine the required precision of the optical align- 
ment, the phase-matching angles were calculated for the optimization 
of different rotational peaks. The results are presented in Table 6.1. 

Table 6.1 


J 

u) R (cm *) 

6(°) 

Y(°) 

e(°) 

8 

75.6 

5.15 

5.11 

0.46 

9 

83.6 

5.42 

5.67 

5.92 

10 

91.5 

5.67 

5.62 

0.56 

11 

99.5 

5.92 

5.85 

0.61 

12 

107.4 

6.15 

6.08 

0.66 


The fact that more than ten rotational lines can be seen simultaneously 
in this large-angle configuration implies that the angle 6 can vary by 
>±1°. This result suggests that the most critical parameter in the 
optical alignment is the overlap of the input beams, and that slight 
errors in the angles are permissible. 
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EXPERIMENTAL DESIGN AND IMPLEMENTATION OF MULTIPOINT CARS 


Many of the factors involved in executing a multipoint CARS experi- 
ment are identical to those described for the single^point CARS work. 

Only the factors that differ will be discussed in this section. 

7.1 Modification of Dye Laser 

The broadband dye laser used for the single-point CARS experiments, 
described in Section 3.1, was also used in the multipoint experiments. 

In order to increase the efficiency of the dye laser, the design was 
modified in several respects. A preamplifier stage was added with the 
pump energy partitioned such that the oscillator was pumped by 10%, the 
preamplifier was pumped by 18%, and the final amplifier was pumped by 
the remaining 72%. The oscillator and the preamplifier stages shared 
the same dye flow system, with the concentration optimized for the 
oscillator. The output coupler was changed from a 4% reflecting wedge 
to a 10% reflecting wedge and the rear mirror was coated to have maximum 
reflectivity at the dye laser wavelength. In addition, all of the dye 
laser optics were AR-coated. The dye cells were changed from 2 cm 
pathlength glass cells to 1 cm fused silica cells in order to optimize 
pumping efficiency. The increase in dye laser energy was approximately 
40%, corresponding to 14 mJ/pulse with 110 mJ of third-harmonic pump 
energy. 

7.2 Optical Arrangement and Alignment 

The optical arrangement for the multipoint CARS is shown in Fig. 7.1. 
In order to achieve the correct angular alignment of the input beams, 
the relative spacing necessary for a given focal length lens was deter- 
mined. Then the appropriate targets and masks were designed and 
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fabricated from aluminum. The parallel input beams were aligned with 
the masks and targets, and then the lenses were inserted. The dye laser 
beam and one of the pump beams were focused with the same 30 cm spheri- 
cal lens. The other pump beam was focused with a spherical lens during 
the initial alignment with the masks and targets. This lens was 
replaced with a 20 cm cylindrical lens for the final alignment. It was 
necessary to have an optical delay in one leg of the pump beam in order 
to have the three beams arrive simultaneously at the interaction region. 
The path of the CARS beam was mimicked by using a HeNe laser, which was 
then imaged onto the slit of the spectrograph. The precision of the 
alignment of the input beams in this manner was not adequate to obtain 
a detectable CARS signal. A fine adjustment to the overlap of the 
beams was performed by using a microscope objective (10X) to view the 
interaction region. After this adjustment at low powers, the laser 
intensity was increased and single-shot burn patterns were made on 
transparent tape in order to confirm complete overlap of the three 
input beams. After this final adjustment, it was possible to obtain a 
detectable CARS signal. In most cases, very little improvement (< 2X) 
in the signal was achieved by tweaking up the detected signal. Most 
often, these adjustments would degrade, rather than increase, the signal. 
This was in sharp contrast to the small-angle single-point geometry, 
where substantial improvements (> 10X) were achieved by optimizing the 
alignment after the signal was initially obtained. It is clear that 
the most critical factor is the overlap of the input beams at their focal 
points, and that slight deviations in the phase-matching angle are 
tolerable. This finding was consistent with the phase-mismatch calcula- 
tions, which also suggested that the overlap was the most critical 
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factor. The optics required realignment each day and, occasionally, 
twice each day. It was often possible to regain the signal by only 
correcting the overlap; however, the signal level tended to be lower 
and eventually it was necessary to realign the entire system. 

7.3 Multichannel Detection System 

The detection system consisted of a 0.5 m spectrograph with vertical 
slits and a 2400 groove/mm holographically ruled grating. Initially, 
the input window of the SIT vidicon was positioned at the exit focal 
plane of the spectrograph. The position was optimized for the maximum 
spectral resolution. Due to the spectrograph astigmatism, this resulted 
in an elongated image which was defocused in the vertical dimension. 

For the single-point CARS work this presented no problems since the 
scanned area of the vidicon could be adjusted to include all of the 
detected signal. However, for the multipoint CARS it was necessary to 
have good resolution in both the horizontal (i.e., spectral resolution) 
and vertical (i.e., spatial resolution) dimensions. The placement of a 
suitable cylindrical lens (30 cm focal length for our system) between 
the final spectrograph mirror and the exit plane brought the image into 
focus for both dimensions. The spatial resolution of this system corre- 
sponded to 0.1 mm in the interaction region. 

Spatially and spectrally resolved multipoint CARS spectra of ^ were 
obtained with this experimental arrangement. However, the signal strength 
was quite low and required the accumulation of a great many laser shots 
to acquire a spectrum which had reasonable signal-to-noise. In order to 
obtain single-shot data it was necessary to increase the magnitude of 
the CARS signal. With the optimization of our NdtYAG laser performance 
and the modifications on the dye laser, we had nearly the maximum input 
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energy that could be expected from our equipment. The optical arrange- 
ment had been optimized with tespect to temporal and spatial overlap 
and energy density in the interaction region. We decided to investigate 
the possibility of amplifying the CARS signal with an image intensifier. 

7.3.1 Image Intensifier 

In order to increase the detected signal level, a 40 mm first- 
generation image intensifier was used in conjunction with the vidicon. 

The image intensifier is a vacuum tube with a self-contained high voltage 
power supply. It consists of an S-20 photocathode, an electron focusing 
arrangement and a P-20 phosphor, which emits green light. The intensi- 
fier is rated as having a white light gain of approximately 80. The 
image intensifier was coupled to the vidicon with a fiber optic rod 
(10 cm x 1.7 cm diam) . This image intensifier /SIT vidicon system, shown 
in Fig. 7.2, will be referred to as the II-SIT detector and the unintensi- 
fied SIT vidicon will be referred to as the SIT detector for the sake of 
brevity. The optical surfaces of the II-SIT system were index matched 
with glycerin, which has an index of refraction close to that of glass 
and which has the advantage of being water soluble and easily removed. 

The alignment of the image intensifier, the fiber optic rod and the 
vidicon is somewhat tedious, since all three pieces must have the coupled 
surfaces parallel to each other. This parallel alignment was accomplished 
by using a HeNe beam reflected from the surfaces. The parallel surfaces 
were brought together by using good translation stages. This method of 
using a fiber optic rod to couple the image intensifier with the vidicon 
was selected because the amount of light loss is lower than would occur 
if lenses were used. The fiber optic coupler was fabricated to be long 
enough so that the vidicon could be used in a temperature cooled mode 
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II-SIT SYSTEM 



Fig. 7.2. Schematic diagram of II-SIT detector system. 
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while the image intensifier remained at room temperature. This is impor- 
tant since the lifetime of the phosphor increases dramatically at lower 
temperatures, causing an unacceptable lag time in the amplified signal. 

7.3.2 Characterization of the SIT and II-SIT Detectors 

Before using the II-SIT, we thought it would be useful to charac- 
terize this detection system and to compare it to the SIT with respect 
to gain and dynamic range. Some preliminary work in our laboratory had 
indicated that these characteristics were different for a pulsed signal 
than for a cw signal, so it was decided that the systems should be 
investigated using the pulsed second-harmonic radiation from the Nd:YAG 
laser. The 532 nm radiation was defocused with a negative spherical 
lens and diffused with a series of ground glass slides. This radiation 
was stopped down with an aperture and focused onto the slit of the 
spectrograph with a 10 cm lens. The intensity of this light was adjusted 
to a reproducible baseline level; this was the maximum light level 
detected. The absolute value of the light intensity was not determined. 
The intensity was variably attenuated with a series of calibrated 
neutral density filters. The illumination of the vidicon was quite 
uniform and a standard scanning format was used to read off the signal. 
The size of the scanned area was 100 x 100 pixels. The single pulse 
signals were reread 20 times in order to insure linearity. The value of 
each different intensity level is the average of 20-50 separate pulses, 
with more pulses being averaged for the lower light intensities. 

The gain achieved with the II-SIT, relative to the SIT, is demon- 
strated in Fig. 7.3. This Is a log-log plot of detected signal level 
as a function of incident laser intensity for both systems. The curves 
cover the complete dynamic range of each detection system. It can be 
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INCIDENT INTENSITY (Arb. Units) 


Fig. 7.3. Log-log plot of detector output as a function of incident 
light intensity for the SIT and II-SIT detectors at 20°C. 





seen that the II-SIT is at its maximum incident light level in the linear 
range just as the SIT approaches its minimum detectable light level. In 
the linear regions of the two curves, the gain is ^50. The similar 
slope of the two curves indicates that the SIT and the II-SIT have a 
comparable dynamic range. The characteristic sigmoidal shape of both 
systems is seen clearly in the analogous linear plots in Figs. 7.4 and 
7.5. At the low end of the dynamic range, an increase in the incident 
intensity does not result in a corresponding increase in the signal detec- 
ted. At the high end of the dynamic range the detector shows saturation 
behavior . 

We investigated the possibility of improving our detectability and 
increasing signal-to-noise by cooling the SIT vidicon while keeping the 
image intensifier at room temperature. We repeated the measurements 
for both the SIT and II-SIT at 20°C, O^C, and -20°C. The curves showed 
no striking temperature dependence. We also performed a number of 
measurements of signal-to-noise at a constant incident intensity in the 
linear regimes of the detectors. The standard deviation of the shot-to- 
shot fluctuations was calculated for the SIT and II-SIT at -20°C, -10°C, 
0°C, 10°C, and 20°C. The standard deviation also showed no significant 
temperature dependence. As a result of this findings, we proceeded with 
the multipoint CARS experiment using the uncooled II-SIT. By using the 
II-SIT, along with improvements in the dye laser and optical alignment, 
we were able to generate a multipoint CARS spectrum in a single laser 
pulse. 
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Fig. 7. A. Linear plot of detector output as a function of incident 
light intensity for the SIT detector at 20°C. 
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Fig. 7.5. Linear plot of detector output as a function of incident 
light intensity for the II-SIT detector at 20°C. 




7.4 Normalization Considerations 


The most difficult problem which remains for the instantaneous 
multipoint determination of scalar parameters like temperature and 
concentration is normalization of the signal. This also presented 
difficulty in the single-point CARS measurements, but the problem is 
even more complex for multipoint CARS. 

The multipoint CARS signal is generated with the line oriented 
horizontally (see Fig. 6.2). The slits of the spectrograph are verti- 
cal, requiring the rotation of the CARS signal by 90°. With this 
spectrograph, the scan format is rotated by 90° relative to that used 
in the room and low temperature single-point work. The channels, corre- 
sponding to spectral dispersion, increase (1-499) horizontally from left 
to right. The tracks, corresponding to the spatial dimension, increase 
(1-499) from bottom to top. When the signal is imaged on the vidicon, 
it should appear as an evenly spaced series of vertical lines, with each 
line representing a rotational peak in the spectrum. For the II-SIT 
system, low channel numbers correspond to the red (low J numbers) and 
high channel numbers correspond to the blue (high J numbers) . The verti- 
cal dimension corresponds to different spatial volume elements, or points, 
along the line from which the CARS signal is generated. Ideally, the 
lines should be the same length and should begin and end in the same track 
numbers. Initially, when the CARS signal was imaged in the vidicon, the 
vertical lines began in different tracks, with the entire series slanting 
up toward the blue end. The degree of this distortion was highly depen- 
dent upon the way in which the signal was directed to the slits (e.g., 
slightly moving a turning mirror could change the angle of the slant 
dramatically). This was due to the fact that the CARS signal was already 
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somewhat dispersed (in the vertical dimension) when it was imaged on the 
spectrograph slits. Changing the rotating optic from a Dove prism to a 
beam-steering unit improved this significantly, as did aligning the tilt 
of the lenses exactly perpendicular to the beam axis. In addition, the 
lines were not all the same length, but tended to become shorter from 
one end to the other. This may have been partially due to the fan-shaped 
astigmatism in the spectrograph. These difficulties were adequately 
resolved for the purposes of demonstrating the occurence of different 
species along the interaction line. However, these distortions would 
have to be reduced further to obtain reliable quantitative results. 

In addition, the pincushion distortion in the vidicon affects the 
interpretation of the spectral and spatial data. This distortion causes 
a straight line of light imaged on the vidicon to appear curved. This 
curvature is most severe at the edges, so that a grid imaged on the 
vidicon would have the appearance of a pincushion. Unfortunately, the 
distortion is not symmetric, so that the correct spatial correspondence 
of the input and output has to be determined for each portion of the 
vidicon. For the single-point CARS, this was not a problem because of 
the normalization procedure chosen. For the CARS in a flame, limited 
portions of the vidicon were used so that the sample and nonresonant 
reference could be imaged on adjacent portions that had the least degree 
of distortion. These alternative approaches are not possible for the 
multipoint CARS since the full two-dimensional capabilities of the 
vidicon are required. A complete mapping of the spatial characteristics 
of the vidicon and a method of correcting the data would be necessary to 
achieve accurate information. A two-dimensional photodiode array would 
solve this difficulty, as well as eliminate the problem of signal lag 
discussed in Section 3. 
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We considered the use of a small-angle three-dimensional single- 
point reference with the large-angle multipoint CARS experiment. We 
anticipate two primary difficulties with this approach. First, the 
single-point reference scheme would not compensate for the effect of 
either permanent or fluctuating spatial inhomogeneities in the sheet of 
pump radiation. Second, because the intensity of the generated CARS 
signal depends on the extent to which phase-matching is achieved, the 
differences in coherence length and interaction volume for the two 
configurations could differentially affect the outlying rotational peaks 
for which the mismatch is the greatest. The influence of this mismatch 
is much more significant for the large-angle geometry than for the 
small-angle geometry. A multipoint reference could be used, however, 
space requirements present a problem. The multipoint CARS experiment 
with no reference requires almost all of the available space on a 
4 X 8 ft. optical table. In addition, the number of precision mounts 
and high-power optics needed to include a multipoint reference is 
virtually doubled. 



8. EXPERIMENTAL RESULTS FOR SINGLE-SHOT MULTIPOINT CARS 


The multipoint CARS spectra presented in the following sections were 
obtained using the large-angle three-dimensional phase-matching geometry 
described in Section 6 and the II-SIT/OMA detection system described in 
Section 7. 

8.1 Spatial Resolution 

The size of the interaction region observed in multipoint CARS 
depends on a number of characteristics of the input beams: (1) the beam 

diameters in the interaction region, (2) the extent of overlap in the 
interaction region, (3) the overall beam energies in the interaction 
region, and (4) the degree to which phase-matching has been achieved. 
Since it is difficult to precisely control all of these parameters, the 
observed interaction region varies somewhat. However, a fairly reliable 
estimate of the interaction region can be made for the present system. 

Laser burn patterns indicate that the spherically focused Stokes 
and pump beam diameters are M).l mm and the cylindrically focused pump 
beam diameter is M).l mm x 2 mm. The beam profiles also were measured 
at lower energies by using a 512 element photodiode array. The inter- 
action region for the large-angle three-dimensional geometry is a 
cylindrical volume approximately 0.1 mm in diameter and 2mm in length 
(see Fig. 6.2). This estimate is confirmed by the results obtained from 
the image size of the multipoint CARS signal on the vidicon. The CARS 
signal was demagnified by 5X using the focusing optics and 0.8X by the 
image intensifier for a combined demagnification on the face of the 
vidicon of approximately 4X. Consequently, an image that has a height 
of one 0.025 mm pixel should correspond to 0.1 mm in the interaction 
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region. It was found that the multipoint CARS image was 20 pixels high 
corresponding to an interaction ^2 mm in length. This result is consis- 
tent with the length of the interaction region determined by direct 
measurement of the input beam dimensions. 

The scan format used for the multipoint spectra divided the image on 
the vidicon into 20 separate tracks of data, so that 20 individual spectra 
are contained in each multipoint spectrum. Each track corresponds to one 
pixel or 0.1 mm, so that all 20 tracks represent 2 mm. A multipoint CARS 
spectrum of Ng is shown in Fig. 8.1. As can be seen in the spectrum, the 
CARS signal extends over all of the tracks. The first and last tracks 
of data frequently have an unacceptable level of noise resulting from 
the way in which the signal is read off the vidicon. For this reason, 
most of the multipoint spectra shown have only 17-18 tracks of data 
plotted, although all 20 tracks were recorded. Consequently, at least 
two additional tracks should be scanned and subsequently discarded in 
the analysis of the data. Because of the beam profile of the sheet of 
pump radiation and because of the diminished overlap at the edges, the 
intensity of the CARS signal tends to decrease at the outlying tracks. 

This is not a problem per se for temperature determination since the 
relative distribution of the rotational lines is unchanged by the total 
decrease in intensity. For concentration measurements, the overall 
decrease in intensity has to be taken into account. 

8.2 Single Species Multipoint CARS 

The spectrum shown in Fig. 8.1 was obtained by passing a flow of 
room temperature ^ gas vertically through the interaction region and 
recording the results of a single laser shot. The asymmetric spike in 
the spectrum is caused by electrical interference from the pulsed laser. 
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Fig. 8.1. Single-shot multipoint rotational CARS spectrum of room 
temperature N£. The spatial resolution is 0.1 mm 
(1 pixel/track) . 
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An example of a multlshot spectrum of Nj Is shown in Fig. 8.2, where there 
are a number of spikes caused by multiple laser pulses. This particular 
spectrum had a different scan format than usual in that each track repre- 
sents 2 pixels, and hence M3.2 mm, in the interaction region. This 
spectrum clearly shows the decrease in the signal intensity at the edges 
of the interaction region, which is located mainly in the middle 10 
tracks. As can be seen in comparing the two multipoint spectra, the 
overall signal-to-noise is substantially better in the multishot spectrum 
than in the single-shot spectrum. 

The single-shot multipoint spectrum of O 2 using the original scan 
format (i.e., 1 pixel/track) is presented in Fig. 8.3. The signal-to- 
noise in this spectrum is better than that in the single-shot ^ spectrum 
because the rotational Raman cross section for O 2 is almost 3X greater 
than the cross section for ^ (Penney et al., 1974). A small nozzle was 
used to produce this flow, so the O 2 concentration is greater in the 
center tracks. This is in contrast to Fig. 8.1 where a diffuse flow of 
N 2 was used, so that the region of gas flow was much larger than the 
interaction volume probed. 

8.3 Different Species Multipoint CARS 

Different species were introduced simultaneously into the interaction 
region in order to determine the ability of the multipoint CARS technique 
to spatially resolve different volume elements along a line. This was 
accomplished by using a coflowing jet configuration consisting of a 1 mm 
diam nozzle surrounded by a coaxial honeycomb. This jet was positioned 
so that the nozzle gas was discharged into the center of the interaction 
region. The spectrum of the nozzle gas is expected to occur primarily 
in the central 10 tracks (corresponding to 'V'l mm in the interaction 
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Fig. 8.2. Multishot multipoint rotational CARS spectrum of room 

temperature N 2 . Spatial resolution is 0.2 mm (2 pixels/ 
track) . Asymmetric spikes are caused by multiple laser 
pulses. 
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Fig. 8.3. Single-shot multipoint rotational CARS spectrum of room 
temperature O2. Spatial resolution is 0.1 mm (1 pixel/ 
track) . 
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region) . The spectrum of the coflowing gas should appear in the outermost 
tracks at either end. Several combinations of nozzle and coflowing gases 
were examined: Nj and 0^* which have a rotational spectrum occurring 

in the same frequency range, and Ar and CCl^F^ (Freon-12) which have no 
spectrum in this range. 

Fig. 8.4(a) shows the single-shot spectrum with as the nozzle gas 
and Ar as the coflowing gas (i.e., Ar-^-Ar) . As expected, the spec- 
trum is primarily in the central tracks, indicating that the CARS signal 
is being generated from the 1 mm diam jet of ^ and the spatial integrity 
of the signal is preserved throughout the detection system. The reverse 
configuration, ^-Ar-Nj, with Ar as the nozzle gas is shown in Fig. 8.4(b). 
The central tracks show no spectrum while the intensity of the N 2 spectrum 
gradually increases toward both ends. 

The effect of a significant density gradient on the multipoint CARS 
spectrum was investigated by using the coflowing jet configuration with 
N 2 and freon. The molecular weight of freon is ^4X that of while the 
weights of Ar and ^ are similar. It was thought that the density gradi- 
ent might induce "phase-front" distortions in the input beams, leading 
to a greater phase-mismatch and a consequent reduction in the generated 
CARS signal. The freon-^-freon spectrum and the complementary 
^-freon-^ spectrum are shown in Fig. 8.5. In Figs. 8.4(b) and 8.5(b) 
the topmost tracks correspond to the region of ^ which occurs before 
the input beams have passed through either the Ar jet or the freon jet. 
Although the relative intensities of the top and bottom tracks change 
from shot to shot, there is no consistent pattern of .the lower tracks 
showing less intensity than the upper tracks. Similar results were 
found with Oj-freon-Oj* shown in Fig. 8.6. We conclude that a density 
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CARS spectrum of (a) 

- resolution is 0.1 ram. 
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Fig. 8.5. Single-shot multipoint rotational CARS spectrum of (a) 
freon-^-freon and (b) ^-freon-^. Spatial resolution 
is 0.1 mm. 
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gradient of this magnitude does not adversely affect the phase-matching 
in the multipoint CARS. 

The final set of two species CARS experiments was done with Nj and 
O 2 flowing from adjacent nozzles approximately 1 mm in diameter. The 
N 2 -O 2 spectrum is shown in Fig. 8.7(a). This spectrum used a scan format 
of 2 pixels/track so that the entire 20 tracks (18 are shown) represent 
^4 mm in the sample region. There is no signal in the outer tracks since 
they represent regions outside the interaction volume. The O 2 spectrum 
can be seen in the upper central tracks and the ^ spectrum is in the 
lower central tracks, with mixing clearly occurring for the tracks in 
between. The lower tracks (2-10) are expanded V3X in Fig. 8.7(b), making 
the amplitude of the ^ comparable to that of the O 2 peaks in Fig. 8.7(a). 
This is consistent with the fact that the Raman cross section for O 2 is 
V)X that for 

8.4 Different Temperature Multipoint CARS 

The feasibility of using single-shot CARS for the instantaneous 
temperature determination of several volume elements along a line has 
been one of the goals of this research. The multipoint CARS of ^ at 
two different temperatures was obtained by using the coflowing jet con- 
figuration described in the previous section, with cold ^ as the nozzle 
gas and hot ^ as the coflowing gas. For the nozzle gas, Nj was cooled 
by flowing it through copper tubing immersed in liquid N 2 . For the 
coflowing gas, Nj was heated by flowing it through copper tubing wrapped 
with heating tape. This was not an ideal arrangement since the jet 
apparatus was fabricated from metal, so that it was not possible to 
maintain a large temperature differential between the two gases. The 
visual difference between rotational spectrum is difficult to determine 
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RAMAN SHIFT 


Fig. 8.7. Multishot multipoint rotational CARS spectrum of (a) 
N 2 -0 2 (18 tracks) with N 2 in the lower central tracks 
and 0 2 in the upper central tracks. Spatial resolution 
is 0.2 mm. (b) The lower 9 tracks (mostly N 2 ) expanded 
-V3X. 





for temperatures that differ by less than 100-200 K. This difficulty is 
compounded in an unnormalized spectrum where both systematic and shot-to- 
shot intensity differences in the input beams influence the intensity 
distribution in the rotational lines tremendously. As demonstrated in 
the single-point CARS experiments, it is essential to have a reference 
leg to normalize a single-shot spectrum in order to even estimate a 
temperature. 

With these considerations in mind, we wanted to demonstrate the 
multipoint single-shot CARS of at different temperatures. The hot 
N 2 gas flowing through the nozzle (no cold N 2 ) is shown in Fig. 8.8(b). 

The cold N 2 flowing through the nozzle (no hot N 2 > is shown in Fig. 8.8(a). 
The differences in these unnormalized spectra are not particularly 
striking, but it can be seen that the spectrum of the hot N 2 is somewhat 
shifted to higher J numbers (larger Raman shift) as expected. The 
combination spectrum, with cold N 2 as the nozzle gas and hot N 2 as the 
coflowing gas is shown in Fig. 8.9. Although the spectrum is quite noisy 
and the signal in the outer tracks is relatively weak, it can be seen 
that the spectra for these outer tracks is shifted to higher J numbers 
relative to the spectra for the central tracks, indicating that the 
surrounding N 2 has a higher temperature than the nozzle gas. The relative 
spectral shift can be seen more clearly in Fig. 8.10, where successive 
tracks are aligned vertically. The spectrum of the warmer gas is also 
less intense because of the decrease in both the number density and the 
population difference factor. It is clear that normalization by a 
reference and a more substantial temperature gradient would produce a 
more dramatic multipoint single-shot CARS spectrum. 
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(a) 




J= 4 6 8 10 12 14 16 

RAMAN SHIFT ► 


Fig. 8.8. Single-shot multipoint rotational CARS spectrum of (a) 
cold N 2 and (b) hot N 2 . Spatial resolution is 0.1 mm. 
Note that the intensity distribution in the hot N 2 
spectrum is shifted to higher J numbers (greater Raman 
shift) . 
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J = 4 
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RAMAN SHIFT — ► 


Fig. 8.9. Single-shot multipoint rotational CARS spectrum of 

hot N 2 ~cold N 2 ~hot N 2 * The intensity distribution in 
the outer tracks is slightly shifted to higher J num- 
bers (greater Raman shift). Spatial resolution is 
0.1 mm. 
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Fig. 8.10. Single-shot multipoint rotational CARS spectrum of 

hot N 2 ~cold N 2 ~hot N 2 » This is the same spectrum as 
shown in Fig. 8.9, but plotted with the peaks aligned 
vertically. 
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CONCLUSION 


One of the initial objectives of this research was to investigate the 
use of CARS as a diagnostic technique for the determination of rotational 
temperatures. The results of the single-point experiments, presented in 
Section 5.1, indicate that CARS provides a means of accurately measuring 
the rotational temperature of N^ at room temperature and below. The good 
agreement between the CARS temperature and the thermocouple temperature 
suggests that the theory used to calculate the integrated intensity distri- 
butions is adequate for this application. The normalization of the 
single-shot spectrum was accomplished by using an identical reference at 
known temperature and pressure. However, this normalization scheme is 
applicable only when there is sufficient overlap of the sample and 
reference spectral distributions (i.e., when the temperature difference 
between the two systems is not too great). The most significant errors 
in the single-point measurements were introduced by the SIT vidicon 
detector. Many of these difficulties could be circumvented by the use 
of an improved detector such as an intensified photodiode array. The 
simultaneous determination of temperature and pressure, discussed in 
Section 5.1.7, requires the use of an iterative procedure because of the 
dependence of the signal amplitude on both temperature and pressure. 

Although pure rotational CARS has been recognized as a means of 
temperature determination for temperatures below 1000 K, it has been 
assumed that rotational CARS has much smaller signal strengths than vibra- 
tional CARS at flame temperatures. We determined the relative signal 
strengths from Ng in a flame for the two techniques and found that the 
peak rotational CARS intensity was 1-2X the peak vibrational CARS intensity. 
The primary difficulty with the rotational CARS measurement was the 
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necessity of normalizing by a nonresonant reference. Because the flame 
temperatures are much higher (1600-1900 K) the sample spectrum cannot be 
normalized by an identical reference at room temperature. Nevertheless, 
if the problems associated with normalization are solved, pure rotational 
CARS is a viable technique for measuring temperatures in flames, and 
would be particularly advantageous for systems requiring the measurement 
of a large range of temperatures (77-1900 K) with the same experimental 
apparatus. 

The ultimate objective of this research was to investigate the 
feasibility of using rotational CARS to determine instantaneous tempera- 
ture and concentration from many spatially resolved points along a line 
intersecting a turbulent flow. The initial obstacle was the relatively 
small signal generated in the large-angle, three-dimensional phase-matching 
configuration. The interaction volume was ^30X less than that for the 
small-angle configuration, resulting in proportionately less generated 
signal in the large-angle case. In addition, one of the three input 
beams was focused to a sheet, rather than to a point, leading to a decrease 
in energy density of ^20X. 

The multipoint CARS signal intensity was too weak to permit the acqui- 
sition of a multipoint spectrum in a single laser shot. Modifications in 
the dye laser provided a 40% increase in output energy. Refinements in 
the optical alignment procedure and optimization of the optical arrange- 
ment also led to increased CARS intensity. The most significant improve- 
ment in the experimental design was the addition of an image intensifier 
to the existing detection system. The new detection- system, referred to 
as the II-SIT, consists of an image intensifier tube coupled to the SIT 
vidicon with a specially fabricated fiber optic rod. The response 
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characteristics of the SIT and of the II-SIT were investigated at different 
detector temperatures. For the detection of pulsed signals, there was no 
improvement in signal-to-noise at lower detector temperatures. The overall 
gain in signal intensity with the II-SIT was ^50. 

With the resulting increase in magnitude of the CARS signal it was 
possible to obtain a multipoint spectrum, consisting of 20 spatially- 
resolved volume elements (0.1 mm x 0.1 mm x 0.1 mm), in a single 10 ns 
laser shot. The CARS signal generated from Nj or O 2 under these condi- 
tions is certainly intense enough to permit the determination of scalar 
parameters (temperature and/or concentration) in a turbulent flow. The 
spatial resolution with respect to position along the multipoint line is 
quite good; for the present system it is ^.1 mm and could be improved 
with different imaging optics. The multipoint CARS experiments vividly 
demonstrate the presence or absence of different species at specified 
positions in the interaction volume, as was shown with the coflowing jet 
configuration using various combinations of Ar, and freon. The 

greatest difficulty to overcome is the normalization of the single-shot 
multipoint spectrum. It is not clear that a single-point reference would 
provide accurate normalization of a multipoint sample. Even a multipoint 
reference, which would increase substantially the complexity of the 
experiment, may not compensate adequately for spatial inhomogeneities in 
the input beams. 

The multipoint CARS technique has the potential for investigating 
more than 20 volume elements along a line. A more powerful laser could 
not improve the present experiment substantially since the present energy 
density is close to the breakdown level of the sample gases. However, 
more energetic input beams could be spread over a larger distance to 
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permit the acquisition of a greater number of volume elements. Since the 
two spherically focused beams (the Stokes beam and one of the pump beams) 
cross at a finite angle, at large distances it would become necessary to 
focus one of these with a cylindrical lens. The interaction volume would 
be determined by the intersection of one line and two sheets of radiation. 

It is clear that the multipoint CARS technique can provide instan- 
taneous, spatially-resolved information regarding scalar parameters in 
turbulent flows. At this point this information is qualitative, and at 
best, semiquant itat ive . If the difficulties associated with normaliza- 
tion of the signal and with distortions introduced by the detector can 
be surmounted, the multipoint CARS technique has the potential of 
acquiring single-shot spatially-resolved quantitative information as well. 
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